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ABSTRACT 

We investigate models of self-consistent chemical enrichment of the intergalactic 
medium (IGM) from z = 6.0 — > 1.5, based on hydrodynamic simulations of structure 
formation that explicitly incorporate outflows from star forming galaxies. Our main 
result is that outflow parameterizations derived from observations of local starburst 
galaxies, in particular momentum-driven wind scenarios, provide the best agreement 
with observations of Civ absorption at z ~ 2 — 5. Such models sufficiently enrich the 
high- 2: IGM to produce a global mass density of C IV absorbers that is relatively invari- 
ant from z — 5.5 — > 1.5, in agreement with observations. This occurs despite continual 
IGM enrichment causing an increase in volume-averaged metallicity by ~ x5 — 10 over 
this redshift range, because energy input accompanying the enriching outflows causes 
a drop in the global ionization fraction of C iv. Comparisons to observed C iv column 
density and linewidth distributions and C iv-based pixel optical depth ratios provide 
significant constraints on wind models. Our best-fitting outflow models show mean 
IGM temperatures only slightly above our no-outflow case, metal filling factors of just 
a few percent with volume-weighted metallicities around 10~ 3 at z ~ 3, significant 
amounts of collisionally-ionized C iv absorption, and a metallicity-density relationship 
that rises rapidly at low overdensities and flattens at higher ones. In general, we find 
that outflow speeds must be high enough to enrich the low-density IGM at early times 
but low enough not to overheat it, and concurrently must significantly suppress early 
star formation while still producing enough early metals. It is therefore non-trivial 
that locally-calibrated momentum-driven wind scenarios naturally yield the desired 
strength and evolution of outflows, and suggest that such models represent a signifi- 
cant step towards understanding the impact of galactic outflows on galaxies and the 
IGM across cosmic time. 

Key words: intergalactic medium, galaxies: formation, galaxies: high-redshift, cos- 
mology: theory, methods: numerical 



1 INTRODUCTION 

The intergalactic mediu m (IGM) contains the majority of 
the Universe's baryons dPave et alJ 1200 lT) and its largest 
structures, yet paradoxically can typically only be painstak- 
ingly observed in absorption against a background source. 
The marriage of high-resolution Echelle spectroscopy on 8- 
10m class telescopes with cosmological hydrodynamic sim- 
ulations in the last decade has led to the understanding 
that the Lyman-a forest arises from highly-photoionized 
Hi tracing fluctuations in the underlying IGM density 
in the so-ca lled Fluctuating Gunn-Peterson Approximation 
JCroft et al This ionization state of the absorbing gas 

is governed by a balance between photo-ionizational heat- 
ing due to the metagalactic UV background and adiabatic 



cooling from Hubble expansion, yielding a tight density- 
temperature rel ation known as the IGM equation of state 
<Hui fc Gnedinlll997D . 

A surprising discovery was that the diffuse IGM is en- 
riched with metals, seen as C iv, Si iv, and O vi absorption 
lines in quasar spectra. Despite our understanding of H I ab- 
sorption, the origin of these metals far from sites of star 
formation remains puzzling. Currently the leading candi- 
date for enriching the IGM is supernova-driven outflows, 
because dynamical disruption of galaxies is too inefficient 
jAguirre et all200ll see also Figure tTUl . Understanding such 
galactic feedback processes is crucial for devel oping a com- 
plete theory of galaxy formation and evolution dDekel fc Silkl 
1986), as observations such as t he shallow faint-e nd slope of 
the galaxy stellar mass function dCole et al .1200 J) versus the 
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halo mass function ( Jenkin^^alJj20Wjh the galaxy mass - 
metallicity relation (iTremonti et al]l2004 lErb et alj f2006'l. 
and the overprodu ction of stars in the Universe in models 
withou t feedback llDave et alJ l200lt ISpringel fc Hernauisd 
l2003bl) indicate that outflows significantly affect galaxy 
properties. 

Observations of metal abundance in various ions that 
trace different densiti es can quantify th e metallicity-density 
relation in the IGM llDave et al.lll998T) . thereby constrain- 
ing the nature of the enriching outflows. Early results from 
Keek's High Resolution Echelle Spectrograph (HIRES) con- 
strained typical IGM metallicities as traced by Civ to b e 
-3 <[C/H]< -2 JSoneaila fc Cowidl99&lDave et alll99St) 
Subse quent observations by ^^^^^Tj^j^I Oscliaveet^ul 
J2003h (hereafter S03), and lSongailaT feoOS) further showed 
that C IV absorption evolves very little from z ~ 5 — > 2, 
which is surprising considering that the vast majority of 
stars in the Universe form at z < 5. These results have 
been interpreted as implying early enrichment by primeval 
galaxies and/or Population III stars at z > 6, where phys- 
ical distances are small and shallow potential wells allow 
winds to distribute meta ls over large c omoving volumes 
([ Scannapieco et all I2002T) . In contrast, lAdelberger et ail 
(2003, 2005) observe enhanced Civ and Ovi absorption in 
the vicinity of galaxies at z ~ 2 — 3, suggesting that en- 
richment is ongoing at lower redshifts in the form of su- 
perwinds from Lyman-break galaxies dPettini et alJ 1200 it 
IShaplev et aill2003f) . Iporciani fc Madaul d2005l) counter this 
by showing that such strong correlations can arise from met- 
als ejected from dwarf galaxies between z « 6 — 12, be- 
cause outflows from highly biased early galaxies should lie 
in highly overdense regions at lower redshifts. Hence the 
sources and epoch of IGM enrichment remain controversial. 

Recent observations of local starbursts have provided 
new insights into the nature and impact of supernova- 
driven winds. Observations of dwarf starbursts and luminous 
infrared galax i es generally find outflows when detectable 
jMartinl l2005l: iRupke. Veilleux fc Sanders! l2005h . observed 
as blueshifted Nal absorption arising from cold clumps en- 
trained in hot metal-rich outflowing gas. While earlier data 
showed litt le trend of outflow properties with host galaxy 
properties dMartinlll999l : Ifleckman et alJl2000h . these more 
recent data suggest that massive galaxies with higher star 
formation rates (SFR's) drive faster and more energetic 
winds. It is of particular interest here that recent data 
exhibit trends that are broadly consistent with theoreti- 
cal exp ectations for momentum-driven (or radiation-driven ) 
winds (iMartinl I200H iMurrav. Quatert. fc Thompsonll2005l) . 
thereby providing for the first time some intuition on the 
physical mechanisms that drive outflows. In such a model, 
radiation from young stars impinges on dust in the outflow, 
which t hen couples to the gas and propels ma tter out of the 
galaxy iMurrav. Quatert. fc T hompson 2005). This is some- 
what different than the canonical scenario where the over- 
pressure from local ISM heatin g causes a bubble t hat even- 
tually bursts out of the galaxy iFuiita et al.ll20o3 . e.g.). An 
advantage of momentum-driven winds is that, unlike heat, 
momentum cannot be radiated away, and hence can plausi- 
bly drive winds over large distances. 

For our purposes, such scenarios for outflows provide a 
physically and observationally-motivated model connecting 
wind properties with host galaxy properties, which we can 



exploit in order to understand the global impact of outflows 
on the IGM across cosmic time. Doing so requires model- 
ing large-scale outflows within the context of hierarchical 
structure formation. This is the main focus of this paper. 

In this paper we explore the impact of various feedback 
mechanisms on IGM metallicity using cosmological hydro- 
dynamic simulations, based on an ex plicit implementation 
of supe rwind feedback pioneered by ISpringel fc Hernauisd 
(2003b) (hereafter SH03). Our simulations also include 
metal-line cooling, which significantly affec ts the tempera - 
ture structure of metal-enriched IGM dAguirre et alJl2005h . 
Our goal here is to construct an enrichment model that 
can reproduce both the metal-line observations of the IGM 
as well as the star formation history of the universe. We 
focus primarily on C IV absorption line observations be- 
tween z = 1.5 — 6.0 since this species is the cleanest 
tracer of IGM metallicity observable over this redshift range 
dSongaila fc Cowielll996tlHellsten et alJll99Sft . 

This paper is organized as follows. §2 discusses the sim- 
ulations including the modifications we made to Gadget-2, 
the various feedback models we run, and the generation of 
simulated spectra. §3 gives an overview of the global IGM 
physical properties in our suite of outflow simulations, while 
§4 discusses the physical characteristics of IGM metals and 
C IV absorbers. §5 compares our findings to observations an- 
alyzed by Voigt profile fitting of lines and the Pixel Optical 
Depth (POD) method. In §6 we show that our results are ro- 
bust against the effects of numerical resolution. We present 
our conclusions in 57. 



2 SIMULATIONS 

2.1 Hydrodynamic Simulations of Structure 
Formation 

W e employ the N-body+hydrodynamic code Gadget- 
2 dSpringell l200Sh . which uses a tree-particle-mesh data 
structure to efficiently compute gravitational forces on a 
system of particles, along with an entropy-conservative 
formulation of Smoothed Particle Hydrodynamics (SPH; 
ISpringel fc Hernauistte'ool to model the pressure forces and 
shocks in the gaseous component. The code is fully adaptive 
in space and time, enabling simulations with large dynamic 
range crucial for studying galaxies together with large-scale 
structure. 

Additionally, Gadget-2 models a number of physical 
processes important for the formation of galaxies. Star for- 
mation is modeled using a subgrid recipe in which each 
gas particle above a critical density where fragmentation 
becomes possible (calculated from the thermal Jeans mass 
based on the local cooling rate) is treated as a set of 
cold clouds embedded in a warm ionized medium, sim- 
ilar to the interstellar medium of our own Galaxy. The 
processes of evaporation and condensation are followed 
analytically within each particle using the formalism of 
iMcKee fc Ostrikerl dl977f) . Stars are allowed to form in the 
cold clouds at a rate proportional to its density squared. 
This produces a disk surface density-star formation rate re- 
lat ionship that is in agreement with the relation observed 
bv lKennicutd 1^998), when a single fr ee parameter, the star 
formation timescale, is set to 2 Gyr dSoringel fc Hernauisd 
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l2003al) . Feedback energy from Type II supernovae is then 
added to the hot phase of the ISM, using an instantaneous 
recycling approximation. SH03 found that the multi-phase 
model produces self-regulated star formation that does not 
suffer from runaway star formation wit h increasing resolu- 
tion, as seen in previous simulations l|Balogh et alJ l200li 
iDave et all 1200 ll) . However, the converged star formation 
rate was still found to be too high when compared with ob- 
servations, motivating SH03 to include galactic outflows as 
we describe in >I2.3I 

Photo-ionization heating is incorporated based on 
a spatially-uniform ultr aviolet background taken from 
lHaardt fc Madaul <200ll) . described in detail in fl2~o1 
Gadget-2 also includes radiative cooling. In the simulations 
run by SH03, the cooling rates were computed assuming pri- 
mordial composition with 76% hydrogen by mass. In >!2.2l we 
describe improvements to the cooling module and extend it 
to include metal-line cooling. 

Gadget-2 tracks the metallicity of gas and stellar par- 
ticles. Gas particles that are eligible to form stars continu- 
ally enrich themselves with metals based on an assumption 
of instantaneous recycling and using a yield of 0.02 (i.e. so- 
lar metallicity). This yield is ro ughly what one would ex - 
pect using supernova yields from lWooslev fc Weaver! (Il995h 
in a Chabrier IMF. When a particle is converted into a star 
(which happens in two stages, such that each star particle 
has half the mass of its original gas particle), then the star 
particle inherits the metallicity of its parent gas particle at 
the time of conversion. 

2.2 Radiative Cooling With Metal Lines 

In the simulations of SH03, cooling is done using an im- 
plicit scheme wherein the cooling timestep is the same as 
the hydrodynamical timestep as set by the Courant condi- 
tion Ispringel fc Hernauisdl2002l) . While in most cases this 
is a stable and accurate method, in dense regions surround- 
ing galaxies the cooling time can be considerably shorter 
than the sound crossing time, but still not so rapid so 
that the particle reaches thermal equilibrium within such 
a time. Since the implicit method uses the cooling rate at 
the end of the timestep to evolve the thermal energy over 
the entire timestep, it can give inaccurate answers in regions 
where the slope of the cooling curve is varying rapidly. This 
happens at moderatel y high densities, and at temperatu res 
~ 10 4 ' 5 - 10 5 ' 5 K fsee lKatz. Weinberg. Hernauisdll996l for 
some examples of cooling curves). 

In order to handle this regime better, we have modified 
Gadget-2 to implement a scheme to follow cooling on the 
cooling timescale. As with lSpringel fc Hernauistl (l2003aTl . we 
assume that the particle's density and non-radiative rate of 
change of thermal energy remain fixed during the dynamical 
timestep; these assumptions are of course not valid in detail, 
but are made to facilitate ease of computation. We also con- 
tinue to assume ionization equilibrium at all times, although 
the algorithm we have implemented makes it easier to fol- 
low non-equilibrium evolution; non-equilibrium effects are 
not expected to be important in the moderately overdense 
IGM that will be the focus of this paper. 

We compute the cooling time as 



where u is the thermal energy, du/dt is the rate of change of 
thermal energy, and e coo i is a tolerance factor that we set to 
0.002. This choice means that in a single cooling timestep 
a particle cannot cool away more than 0.2% of its thermal 
energy. We further limit the cooling timescale such that it 
cannot be less than 0.2% the dynamical timescale, so a par- 
ticle can take up to 500 cooling timesteps for each dynamical 
timestep. The value of this tolerance parameter was chosen 
based on numerous tests on individual particles in the mod- 
erate overdensity, warm-hot temperature regime. 

du/dt for a given particle is obtained from a lookup 
table based on its density and temperature. The lookup ta- 
ble is computed for a given strength of the photo-ionizing 
background, redshift (for Compton cooling), and assum- 
ing primordial composition. The strength of the ionizing 
back ground is int e rpolat ed to the system redshift from the 
lHaardt fc Madaul l200lT) model, and when its strength has 
changed by more than 1% the lookup tables are recom- 
puted. The lookup tables are also recomputed whenever 
Az > 1 since the last lookup table computation, which 
is important in the early universe when Compton cooling 
off microwave background photons is strong. The rate bal- 
ance for primordial species are c alculated as described in 
iKatz. Weinberg. Hernauistl 1^996). If a particle is enriched, 
we add cooling due to metal lines as described below. 

We advance the particle's thermal energy explicitly 
based on the cooling rate computed at the beginning of the 
cooling timestep, and then recompute the cooling rate based 
on the new thermal energy. If the cooling rate has changed 
sign, then the particle has passed an equilibrium point, and 
we return to the original state and reduce the timestep until 
either it no longer changes sign, or the timestep falls be- 
low the minimum value. If the cooling rate has not changed 
sign, we advance the thermal energy using the average of 
the cooling rates calculated at the beginning and end of the 
cooling timestep, thereby preserving second order accuracy. 
We continue to advance the particle until it has been evolved 
over its dynamical (or Courant) timestep. 

To quantify the difference between the old implicit 
method and our new algorithm, consider particles with 
rig « 10~ 3 (i.e. overdensity of 100 at z = 2.5) and T w 
10 4 ' 5 K. By running Gadget-2 with the old and new ver- 
sions of cooling, we found that such a particle has a mean 
absolute difference in the final temperature of 1.4% over a 
single hydrodynamical timestep. At T = 10 5 K it is 0.5%, 
and at higher temperatures it rapidly becomes irrelevant. 
These values are small but non-trivial, and may accumulate 
over many timesteps. At higher densities the effects become 
more significant: For nji = 10 -2 and T = 10 s K, the mean 
difference is 8% (with values discrepant up to ~ x2) and it 
is highly systematic in the sense that the new method pro- 
duces temperatures higher by about 7%. This depends on 
the exact temperature, however, because whether the im- 
plicit method over or underestimates the temperature de- 
pends on the sign of the cooling curve slope at that temper- 
ature. Although the new method increases the total run time 
by typically 20-30%, it seems worthwhile in order to track 
the thermal history of particles more accurately in the mod- 
erately shock-heated, moderately overdense regime, since as 
we shall show (see Figure |5J a substantial amount of C IV 
absorption arises here. 

As Gadget-2 tracks gas-phase metallicities, it is pos- 
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sible to use this information directly to compute the addi- 
tional contribution to cooling rates from metal line cooling. 
To do so, we employ the collis i onal i onization equilibrium 
models of lSutherland fc Dopital lll993t) to generate a lookup 
table of metal cooling versus temperature and metallicity, 
obtained by subtracting their zero- metallicity models from 
their metal-enriched cooling curves. If enriched, a gas parti- 
cle then experiences additional cooling from its metals based 
on a bilinear interpolation within the metal cooling table. 
We also account for the impact of metal cooling on the 
multi-phase subgrid ISM model, since the density at which 
the fragmentation sets in and star formation begins depends 
on the cooling rate. 

2.3 Superwind Feedback 

SH03 found that even with the resolution-converged multi- 
phase ISM model for star formation, the global star forma- 
tion rate predicted by simulations exceeded observations by 
~ x3. Hence they additionally included an explicit model 
for superwind feedback in order to reduce the reservoir of 
gas available for star formation. In Gadget-2, particles that 
are capable of star formation are given a probability of en- 
tering into a superwind based on their current star forma- 
tion rate. If a particle enters a superwind, then it is kicked 
with a velocity given by w w ind, in a direction given by v x a 
(which would yield a polar outflow in the case of a thin 
disk). Furthermore, the wind particle is not allowed to in- 
teract hydrodynamically until it has escaped from the star 
forming region such that its SPH density is less than 10% of 
the critical density for multi-phase collapse; SH03 find that 
the results are not very sensitive to the choice of this value, 
so long as the winds escape the dense star forming region. 
This is intended to mimic a free-flowing chimney of gas ex- 
tending outside the star-forming region as observed in local 
starbursts. 

In the SH03 prescription there are two main free param- 
eters: The wind speed u w ind and the mass loading factor r), 
which is the rate of material being ejected from the galaxy 
relative to its star f ormation rate. Followin g observations by 
iMartinl (Il999l) and lHecknian et alj J2000l) and IGM enrich- 
ment considerations from lAgurrreet^rtJ (120011) . these values 
were both taken to be constant in the runs done by SH03, at 
values of 484 km/s and 2, respectively. We will call this the 
constant wind (cw) model. This model resulted in a z — 
stellar mass density in broad agreement with observations. 

The constant wind prescription, while simple and effec- 
tive, has some deficiencies. Although the nature of the winds 
from the smallest protogalaxies are currently unobtainable 
by observations, using the same large wind velocities for 
these galaxies would heat the IGM too much by z = 3 to 
agree with C IV observations llAguirre et al . 2005), though 
metal-line cooling performed self-consistently may alter that 
conclusion. Additionally, we find that this model has poor 
resolution convergence in terms of global metal enrichment- 
a higher-resolution simulation that resolves small galaxies 
earlier will distribute a great deal more metals throughout 
the IGM at early times as compared to a lower-resolution 
run. 

Observations of outflows from st ar burst g alaxie s have 
improved considerably in recent years. IMartinl J2005T) found 
that the terminal wind velocity scales roughly linearly with 



circular velocity, with top wind s speeds around three times 
the g alaxy's circular velocity. iRupke. Veilleux fc Sanders! 
(2005) studied a large sample of luminous infrared galaxies 
to fi nd that, at leas t when combined with smaller systems 
from IMartinl J^005) , those trends continue to quite large 
systems. It is worth noting that these observations gener- 
ally target cold clouds entrained in the hot wind as traced 
by Nal absorp tion, not the hot wind itself that carries most 
of metals. Still, iHeckmanl J2003I) argues that the wind speeds 
and mass loading factors are likely to be more accurately in- 
ferred from this cold component owing to the observational 
difficulty of detecting X-ray emission from hot gas. 

A feasible physical scenario for the wind driv- 
ing mechanism is deduced by noting that the ob- 
served scaling are wel l expl ained by a momentum-driven 
wind model (IMartinl 2005), such as that outlined in 
iMurrav. Quatert. fc Thompson! (|2005). In such a scenario, 
it is the radiation pressure of the starburst that drives 
the outflow, possibly by transferring momentum to an 
absorptive component (such as dust) that then couples 
to the bulk outflowing material. The presence of large 
amou nts of dust surrounding the classic starburst galaxy 
M82 (jEngelbracht et al.ll200(J) lends circumstantial support 
to this type of scenario. The momentum-driven wind model 
provides us with a physically motivated and observationally 
constrained way to tie outflow properties to the star form- 
ing properties of the host galaxy. Furthermore, it provides 
the impetus for our approach of tying local observations of 
starbursts with outflows across cosmic time. 

The idea of using winds observed in local starbursts as a 
template for winds at all epochs in all galaxies may seem like 
quite a leap of faith. Locally, starbursts are relatively rare 
objects, so it is unclear whether galactic winds are ubiqui- 
tous an d tied solely to the gala xy's star formation rate. How- 
ever, as lHeckman et all ll2000l) showed there is a threshold of 
star formation surface density of around 0.1 M0yr _1 kpc -2 
above which winds are typically seen, and unlike with lo- 
cal star forming disks, virtually all star forming galaxies at 
high redshift satisfy this criterion because they are more 
compact and more vig orously forming stars dHeckman et alJ 
1200(1 lErb et aljl2006r) . Hence down to z ~ 1.5 at least, it is 
plausible that essentially all galaxies are forming stars that 
drive galactic superwinds. Indeed, direct observations o f Ly- 
man break galaxies a t z ~ 2 — 3 bv lPettinT et al] <200 ll) and 
IShaplev et al] ll2003l) show outflows of hundreds of km/s. Of 
course, there is no guarantee that outflows from these high- 
redshift systems follow similar relations as local starbursts, 
but as we shall see this Occam's razor assumption turns out 
to be remarkably successful. 

Guided by local observations, we choose to focus mainly 
on a class of models based on momen tum-driven winds. As 
IMurrav. Quatert. fc Thompson! ll2005l) describes, in such a 
model the wind speed scales as the g alaxy ve l ocity disper- 
sion (see their eqn. 17), as observed bv IMartinl ll2005l) . Since 
in momentum-driven winds the amount of input momentum 
per unit star formation is constant (their eqn. 12), this im- 
plies that the mass loading factor must be inversely propor- 
tional to the velocity dispersion (their eqn. 13). We therefore 
implement the following relations: 

Wwlnd = 3(7 y/ f L - 1, (2) 
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where /l is the luminosity factor, which is the luminos- 
ity of the galaxy in units of the critical (or sometimes 
called Eddington) luminosity of the galaxy, and cro provides 
a normalization for the mass loading factor. Here we as- 
sume that the final radius that the wind is driven to is 
app roximately 100 x the initial rad i us, i.e . r/Ro ~ 100 in 
the iMurrav. Quatert. fc Thompson! i2005T) formalism; the 
dependence on this term is weak. Since we are unable to 
reliably calculate galaxy stellar velocity dispersions directly 
in our simulations owing to a lack of r esolut ion, we approx- 
imate a using virial theorem as a — y — 5$, where $ is the 
gravitational potential at the location of the particle being 
placed into the wind. 

We enlist observations and some theoretical consid- 
erations to determine the fre e para meters cro and /l. 
IMurrav. Quatert. fc Thompson! |(2005) argue that for a 
Salpeter IMF and a typical starburst SED, cro = 300 km/s. 
As we will show, the mass loading factor controls star for- 
mation at early times, so cro can also be set by requiring a 
match to the observed global star formation rate. It turns 
out that for our assumed cosmology, cro = 300 km/s broadly 
matches this constraint as well, so we will use this value 
thr oughout for our m omentum-driven wind models. 

iMartinl d2005Tl sug g ests f L « 2, while 
iRuoke. Veilleux fc Sanders! {2005) find a range of val- 
ues /l w 1.05 — 2. In a momentum-driven wind scenario, 
the relevant luminosity is set by the Lyman continuum 
emission from the stars that is absorbed by the dust 
particles that propel the wind. In such a radiation pressure 
scenario, stars with lower metallicity that produce greater 
Lyman continuum emission (such as those in the early 
universe) would be e x pecte d to drive stronger winds. Stellar 
models by ISchaereil l|2003h suggest an approximate func- 
tional form for far-UV emission as a function of metallicity 
(his equation 1), which we use to obtain the following 
relation: 

f L — j L Q x 1 g-0.0029»(log Z+9) 2 5 +0. 417694 /^v 

For Z — 0.02 (solar metallicity), the value of the exponent 
is zero, making Jl = /l,©, while e.g. for a metallicity of 
10 _3 Zq, Jl = 1.7/z,,0. The metallicity used to determine 
/_l is that of the gas pa rticle en t ering the wind. We typi- 
cally use /l,© = 2 from IMartinl {2005), since the galaxies 
observed in that sample generally have around solar metal- 
licity (or perhaps a tenth-solar, which makes little differ- 
ence) . We will also consider models where we do not include 
this low-metallicity boost, and yet other models where we 
allow the ix,© to randomly vary betwe en 1.05 — 2 in accord 
with iRupke. Veilleux fc Sanded {Hi). 

The form of the mass loading factor 77 is also 
uncertain. While the momentum-dr iven wind model of 
IMurrav. Quatert,. fc Thompson! (|2005l ) predicts 77 oc 1/cr, ob- 
servations seem to find little if any tren d of 77 with a {Martini 
Il999l: IRupke. Veilleux fe Sander j|2o"o5l) . There does appear 
to be a large scatter, so some trends may be hidden in the 
scatter, and it is worth noting that constraining the mass 
loading factor is even more uncertain than determining out- 
flow speeds. Hence we use equation El as our fiducial rela- 
tions, but we will also consider a model where 7J W md varies 
according to equation [5] and r\ is constant. 



Table 1. Simulation parameters 



Name" 


L b 


e c 


m SPH 


mi v 

dark 


M d ' e . 

* ,mm 




w8n256 


8 


0.625 


0.484 


3.15 


15.5 


3.0 


wl6n256 


16 


1.25 


3.87 


25.2 


124 


1.5 


w32n256 


32 


2.5 


31.0 


201 


991 


1.0 



"Additionally, a suffix is added to denote a particular wind 
model as described in i|2.4l 

'Box length of cubic volume, in comoving h~ x Mpc. 
c Equivalent Plummer gravitational softening length, in comoving 
h~ kpc. 

d All masses quoted in units of 10 6 Mq. 
e Minimum resolved galaxy stellar mass. 

Finally, we must circumvent another technical diffi- 
culty with implementing momentum-driven winds within 
a cosmological simulation. The winds are generally driven 
in a manner that continually increases its velocity 
out to some decoupling radius as in equation (15) of 
IMurrav. Quatert. fc Thompson! {2005); this radius can be 
quite large, perhaps ~ 150/i -1 kpc for typical winds. How- 
ever, our implementation only gives an instantaneous kick 
to a gas particle entering a wind, and does not continue to 
accelerate it further. In order to account for this discrepancy, 
we optionally give the winds an additional kick correspond- 
ing to the local escape velocity (given by 2cr), so that the 
escaping wind will have approximately the desired velocity 
as it leaves its galaxy's halo. This may be an overcorrection, 
but if the wind is actively driven to well outside the halo 
scale radius, this should be a fairly good approximation. 

2.4 Runs and Outflow Models 

In summary, we consider the following outflow models: 

• "nw" model: No winds. 

• "cw" model: Constant winds: v w ind = 484 km/s, rj — 
2. This is the model used in SH03 and many subsequent 
papers using those simulations. 

• "zw" model: u w i n d from eqn.|^l a constant 77 = 2, and 
/l from eqn. 4. Note that this is not formally a momentum- 
driven wind model. 

• "mw" model: Momentum-driven winds: w w i n d and 77 
from eqns. 2 & 3, but without a 2cr kick and using Jl = 
/i,0 = 2 (no metallicity dependence). 

• "mzw" model: Momentum-driven winds as above, 
with a varying ft, as given in eqn. 4, /l,q = 2 and with 
a 2cr kick. 

• "vzw" model: Like mzw, but where /l,q is al- 
lowed to randomly vary between 1.05 — 2, as observed by 
IRuoke. Veilleux fc Sanders! d2005l) . 

Table Q lists parameters for our three simulation vol- 
umes having 8, 16, and 32 /i~ 1 Mpc (comoving) box lengths. 
Each volume is run for all of the above wind models. The 
initial conditions used are identical for all the wind models, 
and are generated wh en the universe was still well within the 
linear regime using an lEisenstein fc Hul (Il999t) transfer func- 
tion with f2 m = 0.3, Ha = 0.7, h = 0.7, n = 1, a s = 0.9, and 
Qi = 0.04. Each run has 256 3 gas and dark matter particles, 
with gas particle mass resolutions spanning 5 x 10 Mq to 
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3 x 10 7 Mq. Our smallest volume s till doesn't quite re solve 
the Jeans mass in the high-z IGM JSchave et alJl200(f) . but 
because we will mostly examine C IV absorption which arises 
in moderate overdensity regions, our resolution constraints 
are less stringent, as we discuss in ^ We will mainly focus 
on the 16 /i -1 Mpc box simulations, because as we show in 
our resolution convergence study in SjfjJ this is the largest 
volume for which we can robustly predict C IV observables. 

2.5 Spectral Generation 

From these simulations, we extract the optical depths along 
random lines of sight for 25 ions representing 12 species, 
along with density, temperature, and metallicity for each 
atomic species (which differ because of thermal broadening), 
all as a function of redshift. We use CLOUDY to calculate 
ionization fractions assuming an optically thin slab of gas 
with the given density, temperature, and impinging ionizing 
radiation field (see jg^ . The 25 ions selected have been pre- 
viously observed in optical/UV absorption line spectra, even 
though many ions are too weak to show up in our simulated 
spectra, as well as theoretical expectati ons for IGM densities 
and temperatures jHellsten et al]ll998h . The three ions that 
are considered in this work are H I, C IV, and C ill, while the 
remainder are included only to simulate chance contamina- 
tion. The main contaminants are Nv, Ovi, Sill, Si in, and 
Si iv. When using the pixel optical depth (POD) metho d 
{Ellison et alJll999l: lAguirre et aT1l2002t lAracil et al.ll2004) , 
we must generate the contaminants as accurately as pos- 
sible, and we use alpha-enhanced abundances ([N/C]=-0.7, 
[O/C]=0.5, [Si/C]=0.4 where [C/Fe]=0.0), even though this 
increases the total amount of metals by about 150% because 
half of all metals are oxygen. For much of our analysis how- 
ever, we generate uncontaminated C IV spectra and compare 
to data where contaminants have been removed. Our goal in 
this paper is to explore how simulated observations evolve 
over our redshift range, so we leave for the future a more 
painstakingly detail ed comparison (as done for Lya in e.g. 
iDave fc Trirjpll200ll) . 

Our software, called specexbin, extracts spectra at an- 
gles such that the lines of sight can wrap around the periodic 
simulation box and continuously sample different structure. 
We generally run 30 lines of sight at 30 separate angles be- 
tween 10 and 82 degrees (angles too close to and 90 will 
keep sampling the same structure) for a simulation begin- 
ning at z — 6 and ending at the redshift of the final output 
(usually z = 1.5). We choose 30 lines of sight because we 
want a sample size comparable to a survey of quasar spectra 
achievable currently or in the near-term future, so that we 
can calculate relevant error bars in our plots. A single line of 
sight between z = 6.0 — > 1.5 will traverse a 16 ft _1 Mpc box 
straight across approximately 150 times. We extract from 
simulation snapshots at least every 0.25.Z, while varying the 
ionization background as described in the next subsection 
and smoothly accounting for Hubble expansion. 

specexbin first calculates the physical properties (gas 
density, temperature, metallicity, and velocity) as a function 
of position along the line of sight by averaging the contribu- 
tion of every SPH particle whose smoothing kernel overlaps 
the current pixel. To generate a spectrum in velocity space 
(redshift space), one needs to include the effects of Hubble 
expansion, bulk motions of the gas, and thermal broadening. 



For each ionic species, the ionization fraction is determined 
at each position by using a lookup table where the inputs 
are density and temperature, assuming ionization equilib- 
rium. The ionization fraction is then re-binned into velocity 
space with the inclusion of thermal broadening specific to 
each atomic species. The oscillator strength of the line con- 
verts this value to an optical depth. The metallicity for each 
atomic species in velocity space is saved along with the op- 
tical depths assuming solar metallicity so the we can apply 
any desired metallicity distribution. We typically take the 
metallicity directly from our simulations, but this approach 
retains the option of applying an external p ~ Z relation. 

We then compute fluxes from the optical depths out- 
putted by specexbin. We generate 0.05 A-resolution spec- 
tra, then convolve with an instrumental profile typical of 
an Echelle spectrograph such as HIRES or the Very Large 
Telescope's Ultraviolet and Visual Echelle Spectrograph 
(VLT/UVES), in particular R = X/dX = 43,000, and fi- 
nally add Gaussian noise (S/N = 100 per pixel). As alluded 
to above, we make two type of spectra: C iv-only spectra 
with only the 1548 A C IV component, and complete spec- 
tra with all 25 ions included. We use C iv-only spectra to 
measure C IV column densities and b-parameters, and to 
calculate fi(Clv) (see H5.2.H . We use the complete spec- 
tra when we apply the pixel optical depth (POD) method 
(see H5.3II . because contamination can affect the Civ flux 
decrement. We make multiple spectra for each line of sight, 
placing the quasar redshifts such that all portions of the C IV 
forest from z = 6.0 — > 1.5 are "observed" at wavelengths un- 
contaminated by the Lya forest (e.g. zqso =6.0, 5.0, 4.1, 
3.4, etc.). 

Figure demonstrates how Lya flux, C IV flux, den- 
sity, temperature, and metallicity appear at three different 
redshifts in one of our lines of sight. Some of the trends 
apparent in this model spectrum relate to the main conclu- 
sions of this paper. First, the evolution in C IV absorption 
is significantly less than that seen for Lya due to the com- 
plicated interplay between enrichment, photoionization, and 
IGM heating in expanding large-scale structure. At higher 
redshifts, C IV absorption faithfully traces metal enrichment, 
but at lower redshifts this correlation becomes significantly 
weaker. 

Except where noted, we will group individual Voigt- 
profile fit Civ absorbers (i.e. "components") into "systems" 
if they lie within 100 kms" 1 of another. This is done in 
order to facilitate a more robust comparison among data of 
varying quality and to avoid systematics arising from dif- 
ferent line fitting techniques. For example, in Figure the 
ClV line at z = 1.594 is a single component while the lines 
around z — 4.1 form a multi-component system. 

2.6 UV Background 

We apply a spatially-uniform phot oionizing background is 
taken from lHaardt fc Madaul feOOlf) to the matter distribu- 
tion, both during the simulation runs and during the spectral 
extraction. Unless otherwise noted, the background used is 
the one that is comprised of quasars and 10% of UV pho- 
tons escaping from star-forming galaxies (referred to as QG), 
which turns on at z ~ 9. The slope of the QG background is 
softer than the lHaardt fc Madaul il99rJl quasar-only back- 
ground short-wards of 1050 A (a = —1.8 vs. a = —1.55). 
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Figure 1. Example simulated spectra from the wl6n256vzw simulation at three redshifts, z rs 1.5, 3.0, & 4.5, showing Lya absorption, 
C IV absorption, and the underlying overdensity, temperature and metallicity distributions. Black lines show the properties with peculiar 
velocities and thermal broadening applied, while green lines show the underlying physical distribution. 



The amplitude of the QG background can be con- 
strained by approximately matching the Lya flux decre- 
ment, DLya, to observed values over this redshift range. To 
do so, we divide the background amplitude by 1.6 at all 
redshifts to match the observed mean Lya flux decrement, 
DLya- We do not do this during the simulation run, but 
only in pos t -proc essing as we extract the spectra, but as 
ICroft et alJ <199Sft showed, because photoionization is sub- 
dominant in gas dynamics such a post facto correction yields 
virtually identical results as having done the entire run with 
the lower background. 

In Figure|2]we compar e the mean flux dec r ement in our 
spectra to measurements by Press et al](ll993l) . lRauch et all 
(1997), and iKirkman et alJ (|2005|), whi ch have been cor- 
rected for metal line contamination. The lPress et alJ (I1993T) 
functional fit and Rauch measurements are meant to mea- 
sure the total DLya as they apply corrections for their con- 
tinuum fitting at high redshift, hence we do not continuum- 
fit our spectra for this comparison. Continuum fitting would 
have a non-trivial effect on z > 3 where the Lya trans- 
mission almost never reaches 100% in our spect ra. The dif- 
ferenc e declines at lower redshift; for instance, iDave et alJ 
il997fl find that continuum fitting to the tops of Lya peaks 



results in 5. 7% flux loss at z = 3 an d 1.2% at z = 2. The data 
points from lKirkman et alJ (|2005^> . which have been contin- 
uum fitted, show deviations in the highest redshift bins as 
expected. Their data agree within reason with our measure- 
ments, except at 2 = 3.0 where they find a 22% lower DLya- 
Other than this measurement, we find that the HM01 ioniz- 
ing background simply divided by 1.6 provides a good fit to 
the available data. The fact that we don't need a redshift- 
dependent correction indicates that the HM01 background 
in conjunction with density growth in our assumed ACDM 
universe yields the correct evolution for DLya- 

Interestingly, our outflow models have typically a neg- 
ligible effect on the Lya flux decrements, with the val- 
ues mostly agreeing to 2% among the various simulations. 
This indicates that winds are not affecting the density and 
temperature structure over much of the Universe's volume, 
or put another way, the filling factor of winds is fairly 
small (as we will show more q uantitatively in ^13^31. This 
is consistent with the results of iThe uns et alJ (12002a) and 
iBertone fc White! <2006ft : the latter finding that the only 
quantitative effect of galactic winds on the Lyman-a forest 
is an increase in the the number of sub-A saturated regions 
resulting from dense shells plowed up by the winds. 
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Figure 2. Lyman-alpha flux decrement (Dli/q) averaged over 
30 lines of sight generated from each outflow simulation, using a 
HM01 background amplitude divided by 1.6 across all redshifts. 
This simple correction results in a simulated D/, 9a evolution gen- 
erally in accord with observations. Our outflow models have little 
effect on D^ ya . 



3 GLOBAL PHYSICAL PROPERTIES 

3.1 Star Formation Rate Density 

A major impetus for including superwind feedback is to sup- 
press star formation and solve the overcooling problem. As 
SH03 showed, the constant wind model broadly matches the 
observed star formation history o f the universe, the so-called 
Madau plot iMadau et aljll996l) . Here we examine Madau 
plots for our various wind models to ensure that they fall 
within the observed range as well. 

Figure |3] shows the total star formation rate density as 
a function of redshift averaged over our three simulation vol- 
umes, for all our wind models. D ata points shown are from 
a compilation by iHopkinsl <l2004l) . We include all star for- 
mation in each volume, without selecting out any particular 
galaxy population. In general, all our wind models roughly 
fall within the observed range down to their lowest redshift 
(z = 1.5 or z = 2), with possible discrepancies versus z > 3 
data that are more poorly constrained. With no wind feed- 
back, far too many stars are produced as found by SH03, 
while our various wind models suppress star formation at 
z ~ 3 by ~ x3 — 5 relative to no winds. Hence all our wind 
models broadly pass the Madau plot test. 

More detailed examination reveals that there are some 
differences between wind models that yield insights into 
the impact of varying feedback parameters. For instance, 
in models where we have a variable jj (mzw,mw,vzw) that 
produces high mass loading factors at early times, we end 
up with a global star formation history that is peaked to- 
wards later epochs (z poa k ~ 2 — 4) than in the constant r\ — 2 
cases (cw,zw) which show z pca k ~ 5. This arises because the 
high early mass loading factors in small early galaxies keep 
gas puffy and warm, suppressing early star formation de- 
spite the fact that the wind velocities are low. As a case in 
point, note the vzw and mzw models show virtually identi- 
cal star formation histories at high redshift despite having 
significantly different (mean) values of « W md, showing that 
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Figure 3. The total integra ted star formati on history compared 
to observations compiled by iHo pkins ( 200^) . Wind models sup- 
press star formation from the no-wind case (long-dashed) to bring 
then into broad agreement with data. The mass loading factor 
controls early star formation, while wind speeds become more 
important later on. Outflow models with large rj at early times 
suppress star formation, resulting in a peak in global star forma- 
tion at later times. 

it is indeed r\ that primarily governs early star formation. 
The wind speed is not irrelevant, though, as that represents 
the only difference between vzw and mzw, and they show 
noticeable differences at later epochs. There, the universe 
becomes less dense, winds travel farther, and cooling rates 
are lower, so the wind speed v W md increasingly governs the 
global star formation rate. 

At face value, variable r\ models appear to be in bet- 
ter agreement with observations that show z pca k ~ 2 — 3, 
particularly the vzw model which shows a peak at z ~ 2. 
However, since we are not selecting galaxy populations in 
the simulation analogous to observed samples from which 
the data are calculated, such comparisons are at best pre- 
liminary. For now, we simply note that the feedback recipe 
can have significant impact on the global star formation his- 
tory, including the peak of the star formation rate density 
in the Universe, and all wind models we consider here fall 
broadly within the allowed range. 

3.2 Outflow Properties 

As outflows represent the main new feature of our simula- 
tions, it is worth examining their properties in more detail. 
Here we highlight the differences between wind models in 
their typical speeds, mass loading factors, and energy depo- 
sition rates as a function of redshift. These provide useful 
background information for understanding the behavior of 
C IV absorption in our various outflows scenarios. 

Figure 0] (upper left) shows the mean wind speed as a 
function of redshift for our various wind models. The zw and 
mzw models show similar wind speeds because they both 
employ equation [5] while the mw model shows lowered wind 
speeds primarily because it lacks the 2a kick. The mean 
vzw wind speed is more in line with the mw model at z < 6, 
but vzw produces a wider range of wind velocities due to 
its random assignment of /l,o, and hence ends up being 
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quite different than mw in many other ways. The boost from 
low-metallicity stars (not present in mw) is only important 
at very early times, beca use stars enrich themselves fairly 
quickly (ID aye et all 2006), and once the metallicity exceeds 
around 1% of solar the boost becomes fairly small. The mean 
wind speed when we employ equation|2]becomes comparable 
to the constant wind case (horizontal line) at z ~ 2. It is 
worth noting that the wind speeds at z ~ 3 in all our models 
are comparable to the ra nge of outflow velo cities observed 
in Lyman break galaxies jPettini et aljl200ll) . 

The upper right panel of Figure 2] shows that allowing 
mass loading factor dependence according to equation re- 
sults in much more material being ejected from early galax- 
ies. Of course, because the wind speeds are small and ma- 
terial is infalling gravitationally, the winds do not reach 
very far before being re-accreted. However, they do keep the 
gas somewhat hotter and more diffusely distributed around 
galaxies, so the high mass loading factors are effective at sup - 
pressing early star formation. Note that lDave et alJ (I2006T) 
found that such suppression is necessary in order to match 
observ ed z ~ 6 galaxy luminosity functions. At low red- 
shifts, lErb et all i2006T) estimated a required typical mass 
loading of around 4 in order to understand the trends ob- 
served in the galaxy mass-metallicity relation at z w 2; our 
momentum-driven winds broadly agree with this. 

In the bottom two panels of Figure 0J we show the ki- 
netic energy injected by the winds (left), and the amount 
of energy that reaches the IGM once subtracting off the en- 
ergy required to leave the potential well and enter the IGM 
(right). We normalize these quantities to the average super- 
nova energy for a Salpeter IMF, namely 4 x 10 48 ergs per 
Mq of star formation (the fiducial value used by SH03). 
The cw model was constructed by SH03 to return 100% 
of the supernova energy into kinetic feedback (although 
thermal feedback was included in addition), while the en- 
ergy input from the momentum-driven wind models in- 
creases as galaxies grow. Eventually, the energy inputted 
by winds exceeds the supernova energy; this is physically 
possible because these winds are driven by UV radiation 
from massive stars over their entire lifetimes. An impor- 
tant feature of momentum-driven winds is that deeper 
potential wells in larger galaxies do not inhibit feedback 
of energy into the IGM, in acco rd with observations by 
iRupke. Veilleux fc Sanders! ((2005) . The reverse is true for 
the cw model, where the energy injection into the IGM is 
quite high at high-z and declines to lower redshift, which 
as we sh all see leads to exce ssive hea ting of the IGM (see 
3331 and lAeuirre et alJl2005ft . Indeed. iFerrara et all (l2000l) 
showed that metal enrichment via supernova-driven winds 
alone is energetically insufficient to enrich the IGM to the 
observed levels, and another unknown mechanism was re- 
quired; momentum-driven winds provide such a mechanism. 

3.3 IGM Enrichment and Heating 

Winds add metals and energy simultaneously to the inter- 
galactic medium, and both quantities affect the characteris- 
tics of metal absorption in the diffuse IGM. In this section 
we examine how outflow parameters affect the evolution of 
metallicity and temperature in the IGM. 

Figure|K| upper left, plots the mass-averaged metallicity 
in gas in our various runs. These curves track the cumula- 



tive amount of star formation in each run, with the slight 
difference that metals locked up in stars are not included in 
this plot. Consistent with expectations from Figure |3] the 
no-wind case has the largest amount of metal mass, while 
mzw and vzw have the smallest global metal production. 

Comparing that plot to the volume-averaged metallic- 
ity in the upper right panel shows some interesting differ- 
ences. Note that when quantifying IGM metallicity using 
absorption lines, it most directly traces the volume- averaged 
metallicity. In our outflow models, the amount of metals that 
enters into the IGM depends on a complicated interplay be- 
tween the mass of metals (or stars) formed, the wind speeds, 
and the mass-loading factor. The large wind speeds in the 
constant wind case, together with relatively vigorous star 
formation, produce the most widely distributed metals, with 
a mean metallicity at z ~ 3 of « 1% solar. The difference 
between zw and mzw, which have similar wind speeds, arises 
because the low mass loading factor in zw yields more star 
formation and hence more metals. The lower wind speeds in 
mw and vzw produce relatively low volume- averaged metal- 
licities. Canonical values inferred f rom C IV observations at 
z ~ 3 suggest [C /H1» -2.5 (e.g. ISongaila fc Cowid Il996l: 
iDave et al]ll998l) . which is broadly consistent with all the 
wind models. We will engage in more careful comparisons to 
observations in JS] 

The no-wind model does not distribute metals into the 
IGM hardly at all ([C/H]„ i ume never exceeds -5.0), leav- 
ing most of its volume pristine, and showing that dynam- 
ical stripping of enriched gas owing to interactions can- 
not enrich the diffuse IGM. Our results are in qualita- 
tive agreement with high-resolution Eulerian simulations by 
iGnedin fc Ostrikerl jl997f) with no winds, having spatial res- 
olution comparable to our runs but in much smaller volumes. 
They found a volume- averaged metallicity of«3xl0 _4 at 
z — 4, which is an order of magnitude higher than what we 
find but is still much too low to enrich the diffuse IGM to 
the observed levels. 

The two middle panels show the filling factor with 
metallicity greater than one-thousandth solar (left) and the 
volume-averaged temperature (right). The trends in these 
panels are related. In order to expel metals to large distances 
and fill volume, high wind speeds are required that results in 
high IGM temperatures. Hence the trends are similar among 
our wind models in the two plots: Constant winds produce 
high filling factors and large temperatures, zw and mzw are 
virtually identical, while the mw and vzw models fill less of 
the volume and hardly heat the IGM much above the pri- 
marily photoionization-established temperatures seen in the 
no-wind case. Interestingly, the spread of velocities in the 
vzw model both enrich the IGM while injecting relatively 
little energy, which as we will show turns out to be favored 
by observations. 

The volume filling factor increases with time in all mod- 
els because (1) metals have escaped their potential wells and 
are coasting away from their parent galaxies, and (2) galax- 
ies are forming in less biased regions at later times. The 
volume filled in the cw model begins to asymptote as hier- 
archical buildup occurs and metals ejected by small galax- 
ies falls back into larger halos. The volume will eventually 
asymptote for the other models, but at lower redshift since 
the smaller enriched volumes around galaxies will not over- 
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Figure 4. Wind property averages as a function of redshift for various outflow models. Panels (a) and (b) show the average wind 
velocity and mass loading factor, respectively. Panel (c) shows the energy inputted into the wind particles, and panel (d) shows the 
energy injected into the IGM, after we subtract the energy needed to leave the potential well of the galaxy. The energies are normalized 
by the energy injection from sup ernovae. The wind s p eeds and mass l oading factors at z ~ 2 — 3 in momentum-driven wind models are in 
agreement with observations bv lPettini et al 1 <200j) ; lErb' et alJ l2006fl . The energy input can significantly exceed that from supernovae 
in our momentum-driven wind models, but this is physically plausible since radiation energy is the primary driver of momentum-driven 
winds. 



lap until later, and because the peak star formation is shifted 
to later epochs. 

As a side note. ISchave et ail |2000) demonstrated that 
at z ~ 3, the temperature of the IGM as traced by Lya 
absorption line widths is « 1.5 — 2 x 10 4 K. This tempera- 
ture is above that ex pected for pure photoionization (e.g. 
iHernauist et al.| Il996t or alternatively our no- wind case), 
which they interpret as arising from latent heat owing to 
Hell rei onization. Ou r mzw and zw models broadly agree 
with the lSchave et all <l200(f) measurement, but in our case 
the excess heat is due to outflow energy deposition. Hence 
outflows can produce elevated temperature s without requir- 
ing Hen reionization ('Soka sian et alJl2003T) . We leave more 
detailed studies of Lycv line widths for future work. 

The lower left panel of Figure |S] shows the energy in- 
jected into the IGM per year per /i -1 Mpc 3 , in terms of 
amount of supernova energy produced. This energy injected 
is proportional to 77xSFRxuj GM where uigm is the veloc- 
ity after leaving the potential well of the galaxy. The large 
amounts of early feedback energy in the cw model explains 



how its IGM becomes so hot by z — 6. The momentum wind 
models peak in their energy output at z = 2 — 4, with weaker 
winds peaking later. We will show later that this heat input 
results in significant variations in the global ionization frac- 
tion of C IV, which is a key ingredient in understanding C IV 
evolution. 

The lower right panel shows the amount of metals in the 
IGM where p/p < 100 relative to that at z = 2. This shows 
that in all models, the amount of IGM metals at z — 6 is 
less than one-eighth of that at z = 2. In other words, in all 
of our models the vast majority of metals are injected into 
the IGM during z = 6 ^ 2, rather than at z > 6. Note that 
we do not include any "pre-enrichment" from exotic star 
formation at early times; our models are instead intended 
to test whether normal star formation in ordinary galaxies 
can enrich the IGM with outflow models included. 

Overall, these results highlight the importance of the 
oft-overlooked connection between metal and energy input 
into the IGM. The mass-loading factor r\ roughly governs 
the amount of total metal production at early times, while 
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Figure 5. IGM properties averaged over the simulation box for various outflow models as a function of redshift. Panel (a) displays the 
global mass-weighted gas-phase metallicity, a quantity proportional to stellar mass except that the metals locked up in stars are not 
included. Panels (b) shows the volume-weighted metallicity, which is governed more by outflow speeds than by metal production, as 
the relative trends between models follow that seen in Figure |4ji. Panel (c) shows the fractional volume of the simulation enriched to 
over [Z/H]> —3, which again is broadly governed by outflow speeds. Note that the no-wind cases is unable to enrich more than 0.1% of 
the volume at any epoch. Panel (d) shows the volume-weighted temperature, showing that constant winds dramatically enhances IGM 
temperature, while momentum-driven wind models heat the IGM substantially less. Panel (e) shows the amount of energy reaching the 
IGM in terms of supernova energy per year, per h _1 Mpc; this is essentially a cumulative version of Figure Panel (f) shows the 
amount of metals in the diffuse IGM (p/ p < 100) relative to the z = 2 IGM metallicity; an increase of ~ X 10 from z = 6 — » 2 is typical. 
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the outflow velocity u w i n d mostly determines the physical 
extent to which the diffuse IGM (p/p < 10) is enriched and 
heated. As we shall see, the complicated interplay between 
metal production rate, wind speed, mass loading factor, and 
radiative cooling makes observations of IGM metallicity a 
highly sensitive probe of the physics of large-scale outflows. 



4 PHYSICAL PROPERTIES OF METALS AND 
C IV ABSORBERS 

4.1 Evolution of Metals in the IGM 

Figure |n] shows density, temperature, metallicity, and C IV 
absorption at z =4.5, 3.0, and 1.5, in 15 x 15 ft, - Mpc slices 
that are 100 kms^ 1 wide, from the wl6n256vzw model. 
We have deliberately centered the velocity slice on an over- 
dense region (both in density and C IV absorption) to show 
the large variety of structures formed. This figure illustrates 
some of the trends that will be quantified in upcoming sec- 
tions 1 . 

The growth of large-scale structure, dominated by grav- 
itational instability, is evident in the gas density snapshots 
(top panels) . Outflows increase the metal filling factor while 
growing the metallicity level in previously enriched regions. 
Bubbles of shocked gas grow around star forming systems 
and trace the filamentary structures that house galaxies. In 
the high-redshift snapshot, much of the heating is due to 
winds, as the hot bubbles trace precisely where the metals 
appear. Later, as a proto-group forms in the center of this 
region, the temperature also tracks the virialization of gas 
in th e growing potentia l well. Interestingly, cold mode accre- 
tion jKeres et al ] |2005h is evident along the dense filaments 
feeding galaxies; hence despite the heat input, outflows do 
little to stop rapid accretion at high redshifts, because they 
tend to flow into lower density regions. 

The C IV absorption is shown in front of a backlit screen 
to highlight the morphology of a bsorption, su ch as if our 
Universe was infinite and static j01bersl|l826t) . In reality, 
quasars provide the backlight, and any quasar only probes a 
single pixel. At z = 4.5, nearly all metals show Civ absorp- 
tion when the ionization fractions are highest in the IGM. 
r(Clv)~ 0.05 - 0.10 traces IGM gas with p/p = 3 - 10 and 
t(Civ)> 0.10 traces p/p =10-30 between z = 4.5 - 6.0, 
a time when the vast majority of these overdensities remain 
unenriched by our prescribed winds. By z = 3.0, metals have 
enriched the filaments, which show up as strong absorbers 
connecting the growing galaxy groups. By z = 1.5, much of 
the diffuse IGM carbon is ionized to higher states, and C IV 
preferentially traces higher overdensity structure. 

Figure [7| shows the mass and metal evolution in the dif- 
fuse IGM (p/p < 100, T < 30,000 K) in the top panels, 
the warm-hot IGM (WHIM, p/p < 100, T ^ 30,000 K) in 
the middle, and the condensed IGM (p/p ^ 100, includes 
stars) on the bottom for the cw, mzw, vzw, and nw mod- 
els. The mass fraction in diffuse gas falls in all models from 
2 = 6 — > 2, as mass transitions mostly to the WHIM phase 
resulting from both shock heating on large-scale structure 
and energy feedback. Not surprisingly, the WHIM fraction 

1 For movies of this evolution, see 

http : //luca. as . arizona. edu/~oppen/IGM/ 



follows the mean temperature evolution in Figure The no 
wind model provides a baseline for the amount of WHIM 
formed purely from growth of structure. At z = 2, the cw 
model has a WHIM fraction of 50%, double the nw case, 
indicating half the WHIM results from feedback, and even 
more at higher redshift. The WHIM formed by energy feed- 
back is moderate in the mzw case, and slight in the vzw 
model, thus leading to our conclusion that the fraction of 
WHIM formed by feedback depends most on wind speed. 
Meanwhile, the condensed phase understandably grows the 
fastest in the nw case, while the vzw model clearly distin- 
guishes itself from other two feedback models shown here by 
producing more condensed matter resulting from low wind 
velocities unable to escape their parent haloes. 

Both wind parameters, « w j n d and r/, have a significant 
effect on the fraction of metals in various phases, especially 
compared to the nw case where virtually no metals leave the 
condensed phase. The cw model has a significant amount of 
metals in the condensed phase due to its lower mass loading 
factor, but the metals that do escape into the IGM usually 
are shocked to WHIM temperatures. In momentum-driven 
winds (vzw and mzw) , although rj is larger which suppresses 
star formation, the lower wind speeds allow more metals to 
return to the condensed phase, more so in the vzw case. 

The mean gas metallicities (minus the stars in the con- 
densed phase) are plotted in the right set of Figure |7| The 
condensed gas has a slowly-evolving metallicity of 5-20% 
solar in all wind models at these epochs at z < 4. This 
enrichment level and evolution is similar to that seen for 
damped Lya systems JProchaska et all2003ft . which are ex- 
pected to arise in condensed gas. In contrast, diffuse gas 
shows a steadily increasing metallicity in all models, as star 
formation and winds combine to drive metals into the dif- 
fuse IGM. Interestingly, all wind models show similar diffuse 
phase metallicities and mass fractions, a result of the self- 
regulating nature of feedback (i.e. more feedback curtails 
star formation and metallicity enrichment). The WHIM at 
high redshift is primarily the result of feedback and is more 
enriched than the diffuse IGM. As virialization forms more 
WHIM at lower redshift, the WHIM metallicity becomes 
more similar to that of the diffuse IGM. 

These figures illustrate that winds polluting the IGM 
also affect the temperature structure of the IGM, making 
C IV absorption an evolving tracer of metallicity. Next we 
quantify these trends from the perspective of C IV absorbers. 

4.2 C IV Absorption in Phase Space 

Currently the only observational probe we have of diffuse 
high-redshift IGM gas are one-dimensional quasar absorp- 
tion line spectra skewering a complex matter distribution. 
In this section we present physical properties of the under- 
lying gas for C IV absorption in our simulated spectra, in 
terms of the cosmic phase space (overdensity & tempera- 
ture) of absorbing gas. We focus on the cw, mw, mzw, and 
vzw models as a representative range of wind scenarios: cw 
displays abundant early energy and metal injection, mw has 
comparatively little metals injection and almost no increased 
temperature compared to the no- wind case (see Figure 01, 
while mzw and vzw represent two intermediate cases that, as 
we shall see, are our favored models, with vzw being slightly 
preferred. Hence for illustrative plots we shall utilize the vzw 
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This figure is available in the hill resolution version of the paper at 

http : //nrs a . as . arizona . edu/~r ad /work /carbon/ carbon . pdf . 

Better yet . check out animated versions at 
http : //l nca . as . arizona . ed"u/~oppen/l GM/ . 



Figure 6. A 100 kms" 1 slice of the wl6n256vzw simulation centered on a growing galaxy group at z = 4.5, 3.0, and 1.5. In addition to 
the well- understood growth of large-scale structure (overdensity) , the middle panels show that metal enrichment is tied with the early 
heating of the IGM, especially at high-z before virialization dominates at low redshift. Civ absorption traces nearly all metals at high-z, 
and then declines to trace only the most overdense structures at low-z as the overall ionization correction rises. Close inspection of the 
absorption structures indicates an evolution from smooth and diffuse to clumpy knots tracing the higher overdensities at low-z. 



model, with the understanding that the qualitative trends 
are similar in other models unless otherwise noted. 

Figure |H| shows the ionization fraction as a function 
of density and temperature, in three different redshift in- 
tervals for the wl6n256vzw model. The ionization fraction 
shows two distinct regions: A horizontal band correspond- 
ing to collisional ionization around 10 5 K, and a diagonal 
band at lower temperatures and densities corresponding to 
photoionization. Dashed contours overlaid on top of these 
shaded regions trace the total mass distribution from the 
vzw model. C IV absorption would be visible where these 
contours overlap the shaded regions, if metals are present. 
Metallicity (solid) contours that overlap shaded regions rep- 
resent areas that are in fact traced by C IV for the vzw model. 

One interesting result from this figure is that C IV is 
not an optimal tracer of photo-ionized IGM gas, particu- 
larly at lower redshifts. C V (40. 3A) would be a better tracer, 
but such an observation requires X-ray spectroscopy far be- 
yond the current capabilities. Conversely, there is significant 
amounts of gas heated to around 10 5 K in this and the other 
models, in which C IV should be prominent, even at high 
redshifts. The number of actual observed systems as a func- 
tion of density is governed by an additional factor, namely 
the physical cross-section of the absorbing gas; this will be 
greater at lower densities, which is why weak C IV absorp- 
tion is still more common despite having a relatively low 
C IV ionization fraction. 

Comparing the solid and dashed contours also shows 



that in the diffuse IGM, which holds the majority of baryons 
and the vast majority of volume at these epochs, contains 
few metals. Hence the filling factor of metals is small (as 
seen in Figure|SJ, especially at high-z. The magenta dashed 
line in the two lower redshift bins is the density-temperature 
relationship used by S03 to derive their mass-metallicity re- 
lation, obtained by assuming that the metals occupy the 
diffuse IGM and can be described by a 1-dimensional p — T 
relationship corresponding to the IGM equation of state. 
Their p — T relation is flatter because they assume sub- 
stantial latent heat from H I reionization at z ~ 6; we did 
not include this because latest results from the Wilkinson 
Microwave Anisotropy Probe show that the Universe wa s 
predomi nantly reionized o u t to z ~ 10 iPaee et all 12006^ . 
and our lHaardt fc Madaul i200ll) model assumes reioniza- 
tion at z « 9. The temperatures they use result in a global 
ionization correction that does not change much over their 
redshift range, leading to their conclusion that their lack of 
evolution in the C IV systems implies little redshift evolution 
in the true IGM metallicity. From this plot, it is clear that 
our self-consistent enrichment models predict C IV absorbers 
lying in a wider range of densities and temperature, and are 
not consistent with the assumption that C IV almost always 
arises in a photoionized phase. We quantify this in the next 
section. 

Figure |U] compares the metal mass contours from the 
previous plot (shown as dashed contours) with the distribu- 
tion of the actual C IV absorption (solid) , in the (p, T) plane, 
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Figure 7. The fractional mass, fractional metals, and mass-weighted metallicity of the IGM (from left to right) subdivided into diffuse 
IGM (p/p < 100, T < 30,000 K), warm-hot IGM (WHIM, p/p < 100, T > 30,000 K), and condensed IGM (p/p > 100) (from top to 
bottom). The condensed component includes stellar mass and metals, although the [Z/H] measurement is just for the gas. 



for the cw, mw, and vzw models (top, middle, and bottom 
panels, respectively). We obtain the Civ absorption density 
from the lines of sight in physical space, i.e. without pecu- 
liar velocities and thermal broadening that can significantly 
alter the relationship between absorption and physical prop- 
erties. The two sets of contours illustrate the difference in 
phase space between the true metallicity and the metallicity 
traced by C IV. 

The three wind models shown tell very different stories 
about IGM metallicity, observed and actual. The constant 
wind model has high wind speeds and greater star forma- 
tion at high redshift, resulting in earlier enrichment of lower 
overdensities and the shock-heating of the IGM. The C IV 
ionization fractions grow increasingly smaller at later times 
as the hot metals are ionized to higher states and cannot cool 
efficiently at such low densities. In contrast, the low wind 
speeds of the mw model copiously enrich high overdensities 
while minimizing shock-heating, and lead to a higher overall 
C IV ionization fraction, allowing C IV to more closely trace 
the metals. The intermediate vzw model with its variable 
outflow velocity has both higher wind speeds capable of in- 
jecting metals into the diffuse IGM and lower wind speeds 
replenishing metals in the galactic halo gas. The result is 
IGM enrichment at a wider range of densities than mw with 



less feedback heating than cw, which as we shall see is a 
better match to observations. All forms of feedback result 
in metals occupying a large range of temperatures by z = 3 
resulting in many meta l lines being c ollisionally ionized, a 
prediction also made bv lTheuns et al.l (l2002bl) in their sim- 
ulations with feedback. 

Generally, at high redshifts C IV is a nice tracer of dif- 
fuse metals, but it misses the low-temperature high-density 
regime of the metal distribution (i.e. condensed phase met- 
als) where much of the metals are residing. At z < 3, C IV 
traces diffuse metals less well, because an increasing amount 
of metals arises in hot gas that is too highly collisionally ion- 
ized for C IV. It also fails to trace very diffuse gas that is too 
highly photoionized for C IV. Hence in what could be termed 
the "Age of Civ", this ion between z = 3 — 4.5 is able to 
trace metals over the largest range of overdensities corre- 
sponding to the largest variety of structures. As a side note, 
the relative invariance in C IV absorption across these red- 
shifts can be seen by the similar areas covered by the C IV 
contours, despite the fact that the total metallicity (the area 
of the metallicity contour) is increasing with redshift; we will 
quantify this effect later. 

In the low redshift bin, C IV traces typical overden- 
sities of p/p ~ 100 despite the peak in metallicity mov- 
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Figure 8. Solid (colored) contours show metal mass, while dashed contours show total mass plotted in p — T phase space. Shading 
corresponds to the ionization fraction in Civ. These contours have 0.5 dex steps and are normalized across the three redshift bins. 
The shading goes from — » 40% C IV ionization fraction, where the highest fraction occurs in colli sionally ionized regio ns (T ~ 10 s at 
p/p > 10). The magenta dashed line in the lower two redshift panels is the p — T relation used bv lSchave et alJ l2003t) to derive their 
metallicity-density-redshift relationship. Metals occupy a large region of phase space leading to a variety of ionization conditions that 
evolve with redshift. 



ing downwards to p/p ~ 30. This indicates that Civ ab- 
sorption, particularly strong absorption, is arising primar- 
ily in galactic halo gas. Unfortunately for our simulations, 
if the gas in these halos has a structure similar to our 
Milky Way with cold clumps (i. e. high-velocity clouds) and a 
complex multi-phase st ructure <lMo fc Miralda-Escudell996t 
iMaller &"l3 ulfock 2004T). then we cannot hope to resolve the 
detailed density and temper ature structure in o ur simula- 
tions. Indeed, observations bv lSimcoe et alJ (I2006T) of quasar 
HS1700+6416 between z = 1.8 -> 2.7 suggest sub-kpc 
length scales for many of the strongest absorbers. However, 
lensed quasar image pairs do not show significant devia- 
tions until above a few hundred pc in C IV line profiles in- 
dicatin g highly ionized ga s traces structure on at least kpc 
scales iRauch et aljl200ll) . While these observations probe 
the IGM between z = 1.6 — > 3.6, smaller separations are 
probed at higher redshifts, so it is difficult to determine if 
the characteristic size of the absorber does in fact decrease 
at low redshift as our models predict. 

Figure|n|further shows that C IV arises from both photo- 
ionized and collisionally ionized absorbers; we show a line 
above which collisional ionization dominates. Interestingly, 
the fraction of collisionally ionized C IV is highest at the 
earliest times despite a overall cooler IGM. The diffuse IGM 
densities are too high and the background is too low to ionize 
most metals to C IV, while at the same time the collisionally- 
ionized band intersects regions with metals heated almost 
entirely by early feedback. It is because of this that early 
observations of C IV absorption provide the greatest discrim- 
ination between feedback models, as we shall show more ex- 
plicitly later on. 

The fractions of Civ collisionally ionized between z — 
5.5 — 6.0 are 68% and 43% for the cw and vzw models respec- 
tively. These fractions fall precipitously to 29% and 15% by 



z — 4.0, mainly because the photo-ionized C iv component 
rapidly increases due to the changing ionization conditions. 
The amount of collisionally ionized C IV never really drops, 
and grows rapidly below z = 2.5 in the vzw model as metal- 
enriched gas collapses back into halos (especially overden- 
sities p/p > 300) leading to a ~50% collisionally ionized 
fraction at z ~ 2. 

In summary, C IV absorption in cosmic phase space 
clearly demonstrates that there is no one simple C IV ab- 
sorber type relevant to all redshifts and large-scale struc- 
tures we are investigating. Instead, the photo-ionized and 
collisionally ionized absorption probe two different regions 
of the metal- enriched IGM that evolve independently over 
redshift. As the collisionally ionized component in particular 
is quite sensitive to the outflow model, C IV can in principle 
provide interesting constraints on galactic feedback, espe- 
cially at early epochs. 



4.3 Metallicity-Density Relationship 

Since metals are generated in galaxies that lie at highly bi- 
ased peaks of the early matter distribution (e.g. lDave et al. 
2006), it is expected that there will be a positive correla- 
tion between metallicity and local overdensity. Indeed, such 
a relation was determined based on C IV pixel optical depth 
statistics by S03. They inferred that the median metallic- 
ity followed the relation [C/H] = -3.47 + 0.08 {z - 3) + 
0.65 (log(p/p) - 0.5) between z = 1.8 - 4.1, i.e. that IGM 
metallicity shows a moderate density gradient with little red- 
shift evolution. Such a relationship provides a constraint on 
models of metal injection, though as we shall see this S03 
relationship is sensitive to assumptions about the thermal 
state of the gas which are not supported in our more suc- 
cessful wind models. 
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Figure 9. Dashed contours corresponding to metallicity mass (similar to those in Figure and solid contours corresponding to the 
summed C IV optical depth (derived from lines of sight) in p — T phase space for three different redshift from the wl6n256cw, wl6n256mw, 
and wl6n256vzw models. The blue line divides regions dominated by collisional ionization (top) and photo- ionization (bottom). The cw 
model predominantly enriches the low density IGM while causing significant shock-heating, whereas the mw model enriches mostly higher 
overdensities and heating occurs mostly due to virialization at low redshift. The vzw model manages to enrich a variety of overdensities 
while causing little more shock heating than the mw model. 



The top panels of Figure [TOl show the metallicity-density 
relationship for all our wind models at three redshifts. These 
are not log-linear relationships as have often been assumed 
in fits (e.g. S03). Instead, the metallicity-density relations 
are generally fairly flat out to some moderate overdensity, 
and then drop increasingly rapidly to lower overdensities. 
Furthermore, metallicities are continually growing at all 
overdensities in all our simulations as outflows extend to 
lower overdensities and the condensed regions are constantly 
replenished with new metals from galactic superwinds. 

The cw model shows some different behaviors compared 
to the momentum-driven wind cases. The large early wind 
speeds in cw push metals so far out that the < Z/Zq > 
shows almost no gradient between p/p =1 — 10, while var- 
ious forms of the moment um wind models clearly show a 
gradient. For comparison, ICen et all i2005l) find a "metal- 



licity trough" between p/p — 0.1 — 1 in their TIGER simula- 
tions, because their winds propagate from the most massive 
galaxies anisotropically, enriching the underdense IGM to 
higher metallicities than the filaments. The initial velocities 
in their model are 1469 kms -1 , which quickly slow down 
to a few hundred kms -1 after accumulating material. Our 
cw model is probably most similar to theirs, but with lower 
outflow velocities that are not quite sufficient to produce a 
dramatic inversion in the metallicity-density relation. 

The momentum-driven wind models have a stronger 
metallicity-density gradient, because smaller galaxies living 
at lower overdensities enrich smaller volumes with smaller 
velocities. The vzw model has the smoothest gradient result- 
ing from the assumed spread in velocities enriching a vari- 
ety of overdensities. The zw model, despite its significantly 
lower mass loading factor than mzw, enriches all overden- 
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Figure 10. To p panels show the m ean metallicity as a function of overdensity for all our wind models in three redshift bins. For 
comparison, the ISchave et ail l2003l> relation with its log-normal distribution is included at the two lower redshifts. All models show 
steady growth in metallicity at all overdensities as outflows reach lower overdensities while the metals are continuously replenished 
around galaxies. The bottom panels show the mean Civ/C ratio for all IGM gas independent of metallicity, demonstrating how the 
overall ionization fraction decreases and Civ preferentially traces more overdense gas. Civ absorption will be observed where the 
ionization ratio and the metallicity are both significant. The differences between the models result from the effect of feedback on the 
temperature structures of the IGM. The density-temperature relation assumed by S03 yields a steeper ionization fraction as a function 
of density (solid line), resulting in an inferred metallicity-density relation that is shallower than our simulations. 



sities more efficiently, because it has a higher overall star 
formation rate. Wind speed is the biggest determinant of 
metallicity in the diffuse IGM (especially at p/p < 10) while 
mass loading factor determines the metallicity of the con- 
densed phases of the IGM (p/p > 1000); with larger values 
of r) removing more gas from galaxies, curtailing star forma- 
tion, and decreasing enrichment. These difference emphasize 
the interplay of feedback, star formation, and energy depo- 
sition in the IGM. 

Our metallicity-density relationship is generally higher 
than that inferred by S03 at most overdensities The primary 
reason for this discrepancy is the IGM temperature structure 
in our simulations. As we showed in 21 our simulated IGM 
is continually growing hotter at lower redshift making the 
C IV ionization fraction smaller where the metals are, and 
hence requiring higher true metallicities to match a given 
Q(C iv) absorption. To demonstrate this, the bottom panels 
of Figure [THl plot the overall Civ/C ratio (i.e. the inverse 
of the ionization correction) for all IGM gas independent of 



metallicity as a function of overdensity for the various mod- 
els, along with that derived from the density-temperature 
relation assumed by S03 (the dashed line in Figure |HJ. Civ 
should be primarily observed at overdensities where both 
the metallicity and this ratio are large. 

The S03 relation shows a much steeper ionization frac- 
tion growth with density in the mildly overdense regime than 
any of our simulations. This is because they assume pure 
photoionization, whereas our models have significant contri- 
butions from collisionally ionized C IV. This is even true of 
our no-wind model, owing purely to the growth of structure 
around galaxies. This explains why S03 obtains a shallower 
metallicity-density relation in this regime from C IV obser- 
vations, while (as we will show in we match similar data 
with our steeper metallic ity-density gradient. Interestingly, 
IScannapieco et all <l2006t) suggested that a top-hat distribu- 
tion of metals around galaxies is a better fit to the observed 
C IV auto-correlation function, as compared to the S03 fit. 
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Our simulations produce a distribution that at least comes 
closer to this. 

Fieure lTni is also helpful in visualizing some of the trends 
described in H4.2I First, the Civ/C ratio peak shifts from 
overdensities of about 1 to over 100 from z — 5.5 — > 2.0, 
a startling increase due to the decreasing physical densities 
and (more weakly) the increasing ionizing background. The 
overdensities of the observed C IV do not shift as much be- 
cause lower overdensities have less optical depth, plus the 
diffuse IGM has few or no metals at high-z where the ra- 
tio is highest. Although the photo-ionized and collisionally 
ionized regions are distinct in p — T phase space, the ratio 
merges into a single smooth peak when the temperature axis 
is collapsed (to visualize this look at the shading relative the 
the IGM dashed mass contours in Figure |SJ. Also, the over- 
all decrease in the Civ/C ratio toward low- 2 demonstrates 
the decline in the ability to use C IV as a tracer of metals. 

The differences between outflow models isolate the ef- 
fect of feedback on the temperature structure of the entire 
IGM. At z = 5.5, feedback has not affected much the diffuse 
and underdense gas as can be seen by the similarity in the 
models; only the cw model is able to inject energy signifi- 
cantly to lower the ratio. During the Age of C IV at 2 = 3.5, 
the ionization ratio is high at metal-rich overdensities high 
enough to create significant optical depth in all feedback 
models. The temperature structure causes the most diver- 
gence in the ratio by z = 2 because of the transition to colli- 
sionally ionized C IV, especially in the cw model. Our mod- 
els clearly contrast with the S03 linear density-temperature 
relationship, which intersects higher ionization fractions in 
the photo-ionized regime (see Figure |SJ, leading to their es- 
timates of lower metallicity. 

In short, Civ ionization conditions evolve with den- 
sity and redshift, with an increasing contribution from 
collisionally-ionized absorbers at lower redshifts. Our metal- 
licities at a given density are higher than inferred previously 
by S03 because the broad distribution of metals in the p — T 
plane results in a higher overall C IV ionization correction. 

4.4 Evolution of Median C IV Absorber in Phase 
Space 

Figure fTTI presents the evolution of absorbers in terms of the 
median phase space properties of C IV lines. We show three 
bins of column density (log(AT(C iv)) = 12-13 (weak), 13-14 
(intermediate), 14-15 (strong)) over three different redshift 
ranges (2 = 6^ 4.5, 4.5 — » 3, 3 -—> 2), for three of our wind 
models (cw, mw, vzw) . The values are obtained by summing 
over parcels of gas in redshift space (i.e. including peculiar 
velocities and thermal broadening), and is thereby likely to 
underestimate the density and metallicity of the parcel(s) 
of gas contributing most to the absorption within a given 
line. This is more true at lower redshift as C IV absorbers 
become compressed in comoving space while the comoving 
pathlength per redshift increases. Nevertheless, these out- 
flow models are clearly distinguished, allowing us to draw 
some simple conclusions about the physics of the absorbers. 

In the cw model, the median absorber exists at a tem- 
perature that is dominated by collisional ionization (ss 10 5 
K), which as we shall show produces poor agreement with 
observations. The typical overdensity increases with time 
at all column densities, as diffuse gas becomes increasingly 



ionized to Cv and beyond. In contrast, the mw model's 
absorbers are almost always photo-ionized except at the 
strongest lines at high-2, and show relatively little evolu- 
tion in overdensity over our redshift range, because the low 
wind speeds do not eject metals much outside the relatively 
dense regions that are well-traced by photo-ionized C IV. 
Meanwhile, the vzw model manages to shock heat more met- 
als creating a larger collisionally ionized component, while 
concurrently allowing significant metal enrichment of denser 
halo gas by virtue of some outflow velocities being low. 
This creates a wider range of overdensities and temperatures 
traced by different column densities at any given redshift as 
compared with either the cw or mw models, with weak ab- 
sorbers tracing p/p ~ 4 — 5, and strong absorbers tracing 
p/p « 30 - 100. 

In brief, the behavior of C IV absorption in cosmic phase 
space shows noticeably different behavior for various wind 
models, suggesting that C IV observations can provide a sen- 
sitive probe of galactic outflows. In the next section, we 
demonstrate this explicitly by comparing our outflow models 
to observations. 



5 TESTING OUTFLOW MODELS AGAINST 
OBSERVATIONS 

5.1 C IV Mass Density 

So far we have seen that outflows can have a significant 
impact on the physical properties of the IGM. We now com- 
pare C IV absorption line properties from our various outflow 
models in order to determine which ones, if any, successfully 
reproduce available observations. We begin by examining 
the most global statistic that can be derived from Civ ab- 
sorption, namely the total mass density in C IV absorbers, 
fi(Clv). 

Q,(C iv) provides a benchmark observational test for 
models of IGM enrichment. The observed lack of evolution 
of SI (Civ) from z « 6 — ■> 2 has been cited as evidence for 
very early (2 > 6) enrichment of the IGM |S^ngaiI3|200l|; 
[Scannapieco et all 1 2002; Porci ani fc Ma dau 2005: ISongailal 
2005) from a putative generation of early massive stars 
and/or primeval galaxies. Such a scenario is attractive at 
face value because winds can easily escape the small poten- 
tial wells of early galaxies. But it should be remembered 
that the metals must be expelled to densities approach- 
ing the cosmic mean in order not to have been gravita- 
tionally re-accreted ont o galaxies by lower redshifts (e.g. 
IPorciani fc MadadEjoift . This means winds must not only 
escape galaxies, but must also overcome Hubble flow which 
is quite rapid at these early times, approaching 1000 km/s. A 
population of early galaxies that is energetically able to ex- 
pel winds at ~ 1000 km/s may be difficult to reconcile with 
limits on pre-reionization massive star formatio n from the 
latest WMAP Thompson optical depth results llPaee et alJ 
12006ft . 

None of our models directly test this early IGM enrich- 
ment scenario. Instead, our outflow models attempt to repro- 
duce the observed C iv evolution self-consistently from stan- 
dard stellar populations combined with locally-calibrated 
outflow models. We find that our models are successful, 
thereby removing the need for widespread early IGM en- 
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Figure 11. Median values of Civ absorbers subdivided into column density bins (point size) and redshift (color and shape) along with 
the median values of metallicity ([Z/H], written below with larger text for higher column density). Left, middle and right panels show 
the wl6n256cw, wl6n256mw, and wl6n256vzw models. The medians plotted are derived from spectra that include peculiar velocities 
and thermal broadening. The cw model is dominated by collisionally-ionized absorbers at all redshifts, while mw is photo-ionization 
dominated. The vzw model probes the widest range of densities from weak to strong systems, showing its ability to enrich a wide range 
of environments, while still maintaining mostly photoionization-dominatcd temperatures. 



richment. We further find that observations non-trivially 
con strain our outflo w models. 

ISongailal (1200 lF) determined f2(C iv) by integrating the 
total column density of systems between 10 12 iV(Clv)< 
10 15 cm -2 and dividing by the pathlength using 



n(CIV)(.z) = 



H mc iv £iV tot (CIV,2) 



Cpcrit 



AX(z) 



where 



AX(z) = 



(1 + zfdz 



(5) 



(6) 



These limits are chosen because larger column density sys- 
tems are likely to arise from galactic halo gas, while the lower 
column densities suffer from observational incompleteness. 
When needed, we convert observations to our simulations' 
ACDM cosmology (Qm = 0.3, Qa = 0.7). Measurements 
of Q(C iv) have highlighted some interesting trends that al- 
ready provide a key test for our various feedback models. 

Figure Il2l shows the evolution of f2(dv), as well as 
fi(C), calculated similarly to equation [K] except using the 
total metallicity alo ng our simulated li nes of sight. Observa - 
tions of fl(Civ) bvlSongailal iioOll) . [Pettini et all |2003|) . 
i Bokensberg. Sargent, fc Rauchl (I2003T). ISongailal <2005l) . 
lsimcoeH2006ft . and lRvan- Weber et a!Tll2006l) are s hown We 
also show a measurement at z ~ bv lFrve et alJ (I2003T) for 
reference, although we will not compare to this data here. 

The first striking result from this plot is that while the 
true IGM metallicity continually grows from z = 6 — > 1.5, 
most wind models show fi(Civ) to be relatively constant, 
and perhaps even decreasing at the lowest redshifts. This 
is especially true for our stronger wind models (cw, mzw, 
zw), indicating that outflows play a key role in governing 



C iv evolution. The key point is that in our models, Q(C iv) 
evolution is not a true indicator of metallicity evolution. 

Two counteracting effects conspire to make fi(C iv) ap- 
pear nearly constant. All feedback models grow IGM metal- 
licity at a similar rate, but in the stronger wind models this is 
offset by a decreasing C iv ionization fraction owing mainly 
to an increasing IGM temperature from feedback energy. In 
the case where the feedback energy is minimal such as in 
the mw model, the globally-averaged ionization correction 
(shown in the bottom panel) remains more constant. Nev- 
ertheless, even this model's ionization correction increases 
beyond z = 3 as physical densities of the IGM decrease and 
the ionization background reaches its maximum strength. 
However, for feedback models that heat the IGM signifi- 
cantly (cw,zw,mzw) , the ionization correction from C iv to 
[C/H] increases significantly with time, resulting a fairly flat 
ri(Clv) from z — 5 ^ 2 despite nearly an order of magni- 
tude increase in Q(C). 

Intermediate wind-speed models, mzw and vzw, pro- 
vide the best fit to observations of f2(C iv) over this red- 
shift range. Recall that these are both momentum-driven 
wind models where /l,q = 2 and /l,q = 1.05 — 2, respec- 
tively. Indeed, the vz w model reprodu ces the most recent 
SuperPOD analysis bv lSongailal <l2005l open circles) almost 
exactly, though without a more careful comparison to data 
this should not be taken too seriously. Interestingly, vzw 
and mzw provide good fits to f2(Clv) evolution for some- 
what different reasons. In the mzw model substantial IGM 
enrichment is accompanied by significant heat input, so that 
ionization corrections evolves rapidly to maintain a roughly 
constant $7 (Civ). In the vzw model, there is less IGM en- 
richment owing to lower wind velocities, but also less heat 
input. Indeed, mzw and vzw seem to bracket the allowed 
range of heat /energy input combinations: cw and zw pro- 
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Figure 12. Evolution of f2(dv) (closed circles) and f2(C) (open 
circles), the mass density in Civ and metals, respecti vely. Our 
five fe edback models ar e compared to observ ations from lSongailal 
teOOlfi ( black circles'). IPettini et alJ fcOOaft ( small filled black 
squares ), Ifjoke nsberg. Sargent. &: Rauch| |200 3|) (op en triangles), 
iFrve et alJ fcOOaft ( open square). Isongailal KOOfiB (large open 
dots) ISimcoeH200r3) (large black squares), and lRvan- Weber et alJ 
fcOOrJ) (black triangle is a lower limit). All models show an in- 
creasing f2(C) over this redshift range, but most show a relatively 
constant H(Civ). The ratio is shown as a "global ionization cor- 
rection" in the bottom panel. The cw model is the most metal 
rich, and has increasingly large ionization corrections toward low 
redshift as the IGM is significantly heated and metals are pushed 
predominantly into low density regions. The weak momentum 
wind model (mw) produces little evolution in the ionization cor- 
rection because it does not heat the IGM significantly, and hence 
is the only model that shows a rising Q(C iv) over the entire red- 
shift range. The zw model overproduces f2(C iv) because it pro- 
duces too many metals overall. The mzw and vzw models provide 
good fits to the observations, with the greatest discrimination be- 
tween the two occurring at the highest redshifts. If further obser- 
vations can confirm the lack of a downturn at z > 5, 5, another 
mechanism for early enrichment may be required. 

duce too many metals in the IGM, while mw provides too 
little heat to obtain an ionization correction that grow suf- 
ficiently from 2 = 6^ 1.5. Hence even from a relatively 
blunt tool such as n(Civ), it is possible to place interesting 
constraints on outflow models. 

To further constrain models using fi(C iv) requires 
probing to lower and higher redshifts. Lower redshift C iv 
observations should in principle be able to break this model 
degeneracy, but may be difficult to interpret because C iv 
traces increasingly denser regions that may be subject to 
multi-phase collapse, and the photoionizing background is 



more poorly constrained. Additionally, there is the observa- 
tional challenge of C iv moving into the ultraviolet. Never- 
theless, the dete rmination of fl(C iv) in the local universe by 
IFrve et alJ J2003t) would appear to favor a model with lower 
ionization corrections, since the IGM metallicity should level 
off because of the rapidly dropping global star formation rate 
at z < 1.5. We leave comparisons to z < 1.5 data for future 
work. 

Conversely, moving to higher redshifts may provide 
more robust c onst raints. In particular , new observations by 
ISimcod fcOOd) and lRvan- Weber et alJ fcOOd) at z > 5.5 con- 
tinue to show no apparent drop in fi(C iv) while our models 
predict a gradual decline. These measurements are domi- 
nated by a few strong systems (as these are the only ones 
detectable with current technology), so it is possible that 
may be associated with recent outflows that do not occupy 
a significant volume rather than indicative of widespread en- 
richment. Upcoming near infrared spectra of high- 2 quasars 
discovered by the Sloan Digital Sky Survey should confirm 
just how widespread the C iv is, and whether our models re- 
quire a further enrichment mechanism at high- 2 or perhaps 
an increase in yields from the earlies t stars. The l ack of 
strong C iv observed above z > 5.90 bv lSimco el l|2006 [) alone 
with t he observation of a sudden increase in O i lBecker et alJ 
( 2006) at z > 6 may mean a sudden change in the ionization 
conditions of metals at this epoch corr esponding to the last 
stages of reionization iFan et all ^2006) (R. Simcoe, private 
communication) . 

Finally, we examine a potential source of systematic un- 
certainty in comparing models to observations of f2(Clv), 
namely the effects of Helium reionization. The epoch and 
duration of He n reionization is controver sial, with some ob- 
servations suggesting a rapid transition (lSj3ngailB^_Cowi6 
1996t ISchave et al]l200d: iTheuns et al.ll2002al : IZheng et~ 
2004) at 2 ~ 3, wh ile others favoring a more gradual evo - 
lution before 2 = 3 JPave et alJll998l: ISokasian et alj|2003l) . 
Although the ionization potential of Civ is around 3.5 Ryd 
while Hell is optically thick above 4 Ryd, there may still 
be some residual effect on C iv absorption. To test this, we 
generate spectra using a soft QG ionizing background where 
we have reduced the photoionizing flux above 4 Ry by x 10, 
and derive the resulting 0(C iv) evolution. 

Figure H5I shows the results for the vzw case; the trends 
are similar for other outflow models. C iv traces less over- 
dense material especially between p/p =1 — 10, because of 
the lack of photons capable of ionizing C iv into C V; thus 
decreasing the overall C iv ionization correction. Hence if 
Hen reionization is a sudden process, there should a sud- 
den decrease in fi(C iv) by a factor of around two at the 
epoch. Such a drop is not obvious in the data, although 
error bars are large and could conceivably hide it. At face 
value, it instead appears that He II reionization is a gradual 
process, which could also somewhat affect the evolution of 
O(Clv) possibly making this measure ment e ven more con- 
stant over the redshift range. In fact, 1(1 J2006)l) find evidence 
that the fluctuations in the Hell/Hi ratio require a patchy 
background at energies > 4 Ryd between 2 = 2 — 3. We plan 
to examine the effects of this in more detail in the future. 

To summarize, our models indicate that the invariance 
in fi(C iv) at 2 < z < 5 does not require the majority of IGM 
enrichment to occur at z > 5. Our favored feedback simu- 
lations prefer an interpretation where metals injection fol- 
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Figure 13. The comparison of th e evolution of f2(C iv) with the 
standard lHaardt fc Madaul 1200 ll) ionization background (solid) 
and a softer ionization background where we have reduced the 
ionization at > 4 Ryd by 90% to simulate pre-Hell reionization 
conditions above 2 = 3 (dotted). The same data points described 
in Figure 1121 are displayed. The soft background increases the 
observed C IV by as much as a factor of two. This difference favors 
a gradual He II reionization, as a sudden transition would produce 
an observable drop in Civ. 

lows global stellar mass production, but the accompanying 
energy injection conspires to make 0(C iv) roughly constant 
over the observed redshift range. We cannot rule out an early 
enrichment scenario that invokes a speculative early gener- 
ation of stars to distribute metals in a manner that matches 
SI(Civ) data, but our self-consistent outflow models provide 
a more straightforward explanation. 

5.2 Civ Line Parameters 

5.2.1 Civ Column Density Distributions 

While we showed that fl(C iv) can provide interesting con- 
straints on outflow models, in fact there is more information 
contained in the statistical properties of C iv absorption sys- 
tems from which Sl(Clv) is derived. In this section we com- 
pare our models to observations of the column density distri- 
bution (CDD), which is the number of lines per unit redshift 
interval per log column density bin. The CD D is observed 
to fa i rly accurately follow a steep power law jEllison et alJ 
l200(ilSongaliall200lD . The slope, amplitude, and any devia- 
tions from this power law provide key diagnostics for testing 
models of IGM enrichment. 

We use the AutoVP package foave et al ] l!997ft to fit 
30 simulated spectra per run generated continuously from 
z = 6 — > 1.5. For now we only consider the 1548 A C iv com- 
ponent to avoid confusion; in real spectra, the sparseness of 
C iv absorption generally results in little ambiguity in line 
identification and characterization. This sparseness also pro- 
vides a robust determination of the qu asar continuum , so we 
do not continuum-fit our Civ spectra. iTheuns et al.l (2000) 
demonstrated that AutoVP provided nearly identical fits to 
lines as the standard fitting package VPFIT llCarswell et alJ 
Il987l) for non-saturated lines (as is virtually always the case 
for Civ), and the quality of the fits in the figure appears 



good even in complex systems. One technical difference in 
fitting simulated spectra is that here we fit only the stronger 
component of the doublet, while observations generally re- 
quire the presence of both doublet components to unambigu- 
ously identify C iv. Because of this, observations are more 
likely to miss weak lines. This difference is negligible when 
we integrate the total metal mass, but can become signifi- 
cant when we discuss the small- iV(C iv) end of the CDD. 

To compare our models with observations, we take line 
lists from 9 hig h-resolution spectra taken with H IRES pub- 
lished by iBokensberg. Sargent, fc Rauchl (I2003T) (hereafter 
BSR03). The BSR03 data have typical signal-to-noise ratios 
of S/N ~ 100 and span z — 1.6 — > 4.5. We subdivide the 
BSR03 dataset into systems with z > 3 and z < 3, which 
covers a total of 9.19 and 7.73 in Civ redshift path length, 
respectively. We have also compared with line lists from 19 
VLT/ UVES spectra kindly pr ovided to us by E. Scanna- 
pieco IScannapieco et al]|2006t) . which covers a path length 
of 25.12 from z — 1.5 — 3.0 with lower signal-to-noise than 
BSR03. However, the results are similar to BSR03, and the 
lower S/N and smaller redshift range motivates us to fo- 
cus on the BSR03 data for clarity. In the future we hope 
to obtain the spectra themselves in order to do more careful 
comparisons includi ng detailed noise cha racteristics and res- 
olution effects (e.g. lDave fc TriprJl200ll) . but for now these 
data will suffice in order to ascertain the basic trends versus 
our various wind models. 

Figure 1141 plots CDDs for systems of C iv absorbers in 
three redshift bins, z = 1.5 - 3.0, 3.0 - 4.5, 4.5 - 6.0. We will 
refer to these as the low, intermediate, and high redshift 
bins, respectively, throughout this paper (though we recog- 
nize that z — 1.5 — 3 is not often regarded as "low" redshift). 
For simplicity we omit the zw and nw models from discus- 
sion, because the relevant trends are well-illustrated by the 
other four models (cw, mzw, vzw, and mw). Also shown are 
the BSR03 data (black triangles) in the relevant redshift 
ranges. 

As can be seen, all the simulations generally produce a 
power law between 10 12 ' 5 < N(Cw)< 10 14 ' 5 whose shape 
is in agreement with observations. In fact, at early times 
the simulations show more of a deviation from a power law 
in the sense of fewer weak lines, because the enrichment 
has not spread yet to very low densities; such a trend is 
vaguely noticeable in the data as well. At the very lowest 
column densities the effects of incompleteness become sig- 
nificant. For concreteness, here are our completeness limits 
from simulated spectra with S/N = 100: 75% for 10 12 00 < 
AT(Civ)< 10 12 ' 25 , 88% for 10 12 ' 25 < 2V(Clv)< 10 12 ' 50 , and 
97% for lO 12 ' 50 < iV(Clv)< 10 12 ' 75 . BSR03's completeness 
limits are probably similar but somewhat lower given their 
requirement of identifying the weaker doublet line. As we are 
looking to understand general physical trends rather than 
engage in a detailed comparison versus observations, we do 
not attempt to mimic the data in any greater level of de- 
tail. We simply note that if a simulated CDD exceeds the 
data at the small N(C iv) end, one cannot draw conclusions 
from this without more carefully analyzing the incomplete- 
ness level in the data. 

While qualitatively similar, these four wind models do 
show statistically significant differences, indeed quite large 
ones. The difference in CDD amplitude (and its evolution) 
tells basically the same story as the evolution of fi(Clv): 



22 B. D. Oppenheimer & R. Dave 




Figure 14. Civ Column density distributions for four of our outflow models in three redshift bins, compared to observations by BSR03 
(black triangles). The cw model clearly overpredicts the CDD, while the mw predicts too few weak systems showing that its winds are 
too weak to enrich the low-density IGM. The mzw and vzw model provide promising fits to the data; the overprediction at weak column 
densities is likely a result of lower S/N in the BSR03 data. 



cw produces too much C IV, while mw produces quite little, 
particularly early on. But the decomposition into column 
density bins provides additional information. For instance, 
the mw model does not enrich the low-density IGM much 
owing to its small wind speeds, and hence cannot produce 
enough weak lines at z = 3 — 4.5 to agree with data. We note 
that the incompleteness in the observations is likely greater 
than that in the simulated spectra, so this effect works in the 
wrong direction to reconcile the mw model CDD. So even 
though mw looked promising at early times in the f2(C iv) 
comparison, it actually fails significantly when examining 
CDDs. In contrast, the mzw and vzw models provides very 
good agreement with observations, with the only discrepan- 
cies at iV(Civ)< 10 12 ' 5 which may be explained by larger 
incompleteness in the observations. 

Moving to the low redshift bin, model CDDs draw closer 
for two reasons. First, C IV traces denser gas where met- 
als have been injected earlier on, and second because the 
metallicity in the denser regions starts to saturate, reach- 
ing an equilibrium between galaxies producing new metals 
and outflows driving those metals to lower densities. Again, 
there are still non-trivial differences, as cw still overproduces 
C IV absorption at all column densities, and mw still pro- 
duces fewer weak lines. At face value the mw model seems to 
match the weak lines best, but recall that additional incom- 
pleteness in the data is quite significant here (particularly 
since the BSR03 observations with HIRES are less sensitive 
in the blue), so the turnover in the data at iV(Civ)< 10 13 
may be at least partly artificial. There is a hint that all 
models produce too few of the very strongest systems with 
JV(C iv)> 10 14 ' 5 , which may be due to an increasing popula- 
tion of C IV absorbing gas in small, cold clumps that our sim- 
ulations cannot resolve. For instance, ISimcoe et alJ <l2006t) 
derive small thicknesses (< 1 kpc) and high overdensities 
(p/p > 100) for many of their strongest C IV components, 
which is below our spatial resolution limit. 



Overall, the C IV CDD provides a sensitive test of the 
manner in which outflow models push metals into the IGM. 
Comparisons to observed C IV column density distributions 
favor the vzw and mzw models, which shows little redshift 
evolution at all but the smallest column densities, and whose 
total C IV absorption is close to that observed. More gener- 
ally, it demonstrates that outflows must have wind speeds 
large enough to enrich the very diffuse IGM (giving rise to 
weak C IV absorption) at early times, but must not be so 
strong as to overproduce C IV absorption overall. 



5.2.2 Civ Linewidths 

While column density distributions probe the spatial distri- 
bution C IV in the IGM, line widths are more directly related 
to the temperature and dynamics of the absorbing gas. For 
Hi, the line widt hs at z > 2 are domin ated by Hubble flow 
broadening (e.g. iHernauist et alJll99dh . i.e. the cosmic ex- 
pansion velocity spread across an unbound absorbing struc- 
ture. This becomes less true at lower redshifts, where H I 
line widths obtain an increasing contribution from thermal 
broadening jDave fc Tripdl200lT) . mainly because they arise 
in higher overdensity struc tures that do not have substantial 
Hubble flow across them dPave et al]ll999l) . In the case of 
C IV, even at high redshifts they arise in fairly overdense re- 
gions (see Figure so we can expect that their line widths 
would have a significant contribution from thermal broad- 
ening. 

Figure [TTjI shows histograms of &-parameters (i.e. Voigt 
profile line widths) of components with 10 13 < A r (Clv)< 
10 14 cm -2 in four of our outflow models for our three red- 
shift bins, compared with BSR03 observations (filled circles 
in the lower two redshift bins) . This column density range is 
optimal because observations both yield a statistically sig- 
nificant number of lines and generally approach 100% com- 
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Figure 15. Line width (6-parameter) distributions for 10 13 < iV(C iv)< 10 14 cm -2 systems for our various outflow simulations compared 
to observations by BSR03 (black dots). The arrows at the bottom show to the median values of the distribution for each model, and the 
top arrow is the median of the BSR03 data. The data agree with the mzw model at intermediate redshift, but suggest more compact, 
cooler Civ absorbers at low redshift such as those in the vzw or mw model. 



pleteness. Arrows at the bottom of the plot show the median 
values. 

The trends in this plot follow those in Figure mid- 
dle right panel (which showed the mean temperature evo- 
lution of the IGM). The cw model has significantly larger 
6-parameters owing to substantial shock heat input from 
outflows, and appears to be clearly discrepant when com- 
pared to observations. Indeed, the cw model actually has 
more lines (~ 66%) at z > 4.5 that arise from collisional 
ionized C IV (T > 10 5 K, b therma i > 11.8 kms" 1 ) than 
at 2 < 4.5 where C IV is photoionized in cooler gas; this 
partly explains the b me d evolution downward. Observations 
from BSR03 (solid points) suggest the mzw model is the 
best fit at intermediate redshift, while the vzw model better 
approximates the low-z absorbers. The momentum-driven 
wind models (mw, mzw, vzw) are remarkable in that 6 met j 
does not vary by more than 1 kms -1 from z = 6 ^ 1.5, es- 
pecially considering that one of the main points of this paper 
is that C IV absorbers show significant physical evolution. 

As an informative exercise, we subdivide the width of 
the vzw model line into a thermal component, bth, that 
depends only on temperature, and a spatial component, 
b sp , due to Hubble expansion across the physical extent 
of the absorber. We ignore turbulent broaden ing since the 
IGM at these overdensities is very quiescent iRauch et alJ 
2005). Specifically, we assign b to t = \J b\ h + b 2 sp , where 

b t h — 3.7-\/T/10 4 for carbon, and b sp = H(z) w(Civ) where 
w(C iv) is the spatial extent of the C IV. Figure HrJl shows the 
total (solid) and thermal (dashed) fa-parameter histograms 
for the vzw model, with bth calculated from the tempera- 
tures at the line centers. The thermal widths are the most 
significant component of the total line width, especially at 
low redshift, demonstrating (as expected) that at the mod- 
erate overdensities of C IV absorbers the Hubble broadening 
is sub-dominant. This implies that 6-parameters can effec- 
tively trace the IGM temperature in Civ absorbing gas. 



In summary, the line width distributions and the low- 
column density end of the CDD together put significant con- 
straints on the wind speeds of early outflows: They must be 
strong enough to pollute the diffuse IGM at early times, but 
weak enough so as not to overheat the IGM. Some heating is 
accommodated, and in fact during the Age of C IV there is a 
hint that without wind heating, line widths are too narrow. 
More careful comparisons will be necessary to place more 
detailed constraints, but the fact that our locally-calibrated 
momentum-driven wind scenarios naturally fall within these 
relatively tight constraints represents a significant success 
for these models. 

5.3 Pixel Optical Depth Statistics 

The pixel optical depth (POD) method provides a powerful 
way to measure the metallicity in the diffuse IGM across 
a wide range of densities. The idea is to directly utilize 
optical depth information in the spectrum rather than fit 
lines, so that even weak absorption that would not yield 
a st atistically significant line identification can be quanti- 
fied. ISongaila l <^)98|) pio n eered this approach, which was 
also used bv iDave et al.l l|l99£l) to quantify O VI absorp- 
tion in the crowded high- redshift Lv-q fo rest. More recently, 
lAguirre et alJ (120051) and lSongailal i2005l) have expanded on 
POD analyses to precisely quantify C IV absorption in the 
IGM. In this section we compare POD analyses of simulated 
spectra with such observations. 

In POD analyses, a histograms of optical depth ratios 
of two ions is constructed. This is done by identifying signif- 
icant absorption in one ion, and examining the location in 
the spectrum where the other ion should be, and inferring 
the optical depth ratio of the latter to the former. In our 
case, we will consider C iv/ H I and C iu/C IV rat ios, in order 
to compare with data from lAguirre et all l)2005l) . Since C ill 
has a rest wavelength of 977A, we must also correct for coin- 
cident absorption from H I Lyman series lines; we apply the 
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Figure 16. Histograms of total linewidth (solid lines) for the vzw model and thermal linewidth (6 £ ^) derived from the temperature 
of the Civ absorbing gas (dashed lines). Thermal broadening comprises the majority of the total linewidth, especially at low redshift, 
showing that the fe-parameter distribution can provide a sensitive test of feedback energy input into the IGM. 



method detailed bv lAguirre et all J20021) to do so. Specifi- 
cally, the lowest order unsaturated Lyman-series line is used 
to determine the H I optical depth, which is then subtracted 
from all Lyman lines leaving the absorption component of 
any metal lines. Higher-order Lya lines with their lower os- 
cillator strengths allow one to probe much higher H I, up to 
r(Hl)> 1000 as opposed to t(Hi)~ 3 using Lya alone. The 
POD method is surprisingly effective at extracting metal 
lines in the Lyman forest, although the effectiveness declines 
at z > 4 due to the thick Lya forest. While the full POD 
histograms are rich in information, for brevity here we focus 
on the median of the optical depth ratios as a function of 
Hi or Civ optical depth. 

We apply the POD method to the 30 lines-of-sight 
extracted from each simulation, where we include all 25 
ions and continuum-fit the spectra blueward of Lya as de- 
scrib ed in H2.6I We c ompa re our results with POD statistics 
from lAguirre et alJ i2005tl for a sample of 6 quasars with 
< z >w 3.1. 



5.3.1 Cw/Hi POD Statistics 

Figure H71 (top panels) shows the median Civ/Hi ratio as a 
function of Hi optical depth for a variety of redshift ranges, 
for the cw, mw, mzw, and vzw models (left to right). This 
POD statistic shows even more dramatic differences between 
the various feedback models t han the CDD s. For compar- 
ison, we plot the data from lAguirre et alJ i2005l) in the 
z — 2.8 — 3.4 panels. 

The overall trend is that the median POD ratio rises 
with time. This is primarily because the Hi optical depth 
at a given overdensity is dropping more rapidly with time 
than the C IV optical depth because as the universe becomes 
less dense H I becomes more ionized, while C IV ionizing to 
C v is offset somewhat by C ill ionizing up to C IV. However, 
a secondary non-trivial effect is that the IGM is becoming 
increasingly enriched. While the median POD ratio is rel- 



atively flat with Hi optical depth, the peak does shift to 
lower Hi optical depth with time, owing to winds pushing 
metals further out into the IGM. 

Comparing to observations, it is clear that the cw model 
overpredicts the median POD ratio at z ~ 3.1, consistent 
with over-enrichment of the IGM as has been seen from the 
fi(Clv) and CDD comparisons. What is new here is that the 
shape of the curve is also discrepant, in that the peak value 
is produced around logr(Hl)~ 1, while observations prefer 
a peak at logr(Hl)~ 1.5. This is indicating that metals are 
too widely distributed in the cw model, in addition to being 
overproduced. 

The mw model represents the other extreme, where it 
fits only the high density regions (log r(H l)> 2) but does not 
distribute the metals widely enough to match at logr(Hl)< 
1.5. Clearly, the winds are too weak in this model to enrich 
the diffuse IGM. 

In contrast, the mzw does and excellent job fitting the 
full range of H I optical depths, indicating the winds are both 
producing the correct amount of C IV as well as distributing 
it correctly versus the underlying Hi absorption. The vzw 
model is nearly as good, although it appears to produce 
slightly too little C IV absorption overall. Without doing a 
more careful comparison, we can only say that both these 
mo dels are in reasonable agreement with observations. 

lAguirre et alJ i2005l) applied the POD method on sim- 
ulations including a 10 ft _1 Mpc, 2 x 216 3 particle Gadget- 
2 simulation from SH03 (the Q4 simulation) with the cw 
feedback model, but without metal-line cooling. They added 
an ad hoc cooling method wherein SPH particles with cool- 
ing times less than Hubble times instantaneously reach 
T ~ 10 4 K; resulting in a C iv/H I ratio increasing by about 
~ xlO with cooling. Using a test simulation, an 8 /i _1 Mpc 
box vzw simulation with 2 x 128 3 particles with and without 
metal-line cooling, we find the metal-line cooling increases 
C IV absorption only by between 60-90%, likely because in 
the denser regions the IGM is continuously reheated by feed- 
back, while in the underdense regions the cooling time even 
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Figure 17. Median pixel optical depth ratios for our outflow models in three redshift bins, compared to observations bv lAguirre et all 
12005) (solid lines). The mzw model provides the best fit to the observed Civ/Hi ratio where Hi is a proxy for overdensity, indicating 
that the metals are properly spatially distributed in this model. The vzw model does nearly as well, while the cw model injects far too 
many metals and the mw model is only able to enrich the most overdense regions. The vzw and mw models adequately fit the Cm/Civ 
ratio while the cw and mzw models underestimate this ratio due to excessive heating of the IGM. The mismatch at t(Cih)~ —0.1 is 
likely a result of our simulated spectra having lower S/N than the observations resulting in worse Lyman series subtraction at high 
optical depths in the Cm forest. The Civ/Hi ratio grows toward lower redshift due to increasing metal enrichment, while the Cm/Civ 
decreases due to heating of the IGM and decreasing physical densities. 



with metals approaches the Hubble time. Hence the addi- 
tional effects of metal line cooling cannot reconcile the cw 
model with observations. 



5.3.2 Cm/Civ Statistics 

The Cm/Civ POD relation is used by S03 to provide a 
constraint on the temperature of the absorbing Civ. This 
is a difficult measurement to obtain since Cm (977A) is 
buried in the Ly/3 forest, but their method applied identi- 
cally to observations and simulations shows that the ratio 
can still be instructive. They conclude that the high ratios 
(t(C iii)/t(C iv)> 3 at z > 3) rule out large fractions of 
C iv-absorbing gas at T > 10 5 K where collisional ionization 
dominates and the ratio drops precipitously (refer to their 
Figure 7 to see how this ratio depends on density and tem- 



perature using the QG background). lAguirre et all l l200Eh 
find that their Gadget-2 simulation remained much too 
hot to reproduce the observed ratio, and ratios were 0.65 
dex too low even when including metal cooling, only a mod- 
est improvement. 

We find that our cw model indeed produces too low a 
median POD ratio. This indicates that the IGM is too hot 
in this model compared with observations, as we also saw 
in the 6-parameter comparison. In contrast, our momentum- 
wind models fare considerably better. The mzw model still 
appears to heat the IGM too much, but the vzw and mw 
models (with relatively low wind speeds that do not add 
much heat to the IGM) fare quite well. 

S03 points out that the observed C m/C IV ratios are in- 
consistent with the absorbing gas being at T > 10 5 K. How- 
ever, even the vzw and mw models have a substantial frac- 
tion of their C iv-absorbing gas collisionally ionized, about 
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30% at z ~ 3. Despite this, these models match the POD ra- 
tios quite well, showing that some collisionally-ionized ClV 
is still accommodated by observations. 

S03 further observes a gradual decline in 
t(C iii)/t(C iv) with time. All wind models produce a 
gradual decline qualitatively consistent with these data. 
In photo-ionized gas this decline is mainly driven by 
Hubble expansion increasing the IGM ionization level, and 
secondarily by the increasing temperature of the IGM. 

In summary, pixel optical depth techniques provide a 
complementary technique to Voigt profile fitting for exam- 
ining IGM C IV absorption. The overall results from both 
of these techniques support the same basic conclusions: Our 
cw model overproduces metals, distributes them too widely, 
and makes the IGM too hot, while the mw model under- 
produces metals, does not distribute them widely enough, 
and keeps the IGM (possibly) too cool. In between lie the 
mzw and vzw models that seem to span the allowed range 
of metal and energy input. More generally, these observa- 
tions already provide a non-trivial bracketing of the allowed 
energy and metal input into the IGM as a function of red- 
shift. More careful comparisons with observations, particu- 
larly pushing to higher and lower redshifts, should allow us 
to provide strong constraints on the properties of galactic 
outflows across cosmic time. 



6 NUMERICAL RESOLUTION 

Numerical simulations by necessity have a finite dynamic 
range, and even advanced ones such as ours cannot fully 
resolve the relevant scales from star formation to the Hub- 
ble volume. Hence we must ensure that our results are not 
sensitive to issues of numerical resolution. For this purpose 
we have run with 8, 16 and 32 /i~ 1 Mpc box lengths, as de- 
scribed in which span a factor of 4 in spatial resolution 
and a factor of 64 in mass resolution . The results presented 
so far have been derived from the 16/i _1 Mpc volumes for 
each outflow model, because as we will show in this section, 
their results show good convergence relative to the higher- 
resolution 8/i _1 Mpc runs (at least for the momentum-driven 
wind models), but the 32/i _1 Mpc runs do not show good 
convergence. We present simple arguments that explain this 
behavior in terms of the Jeans length in the moderately- 
overdense IGM giving rise to C IV absorption. Finally, we 
show that the constant wind model shows poor convergence 
at all resolutions. We focus on the convergence properties 
of the observable quantities that we have compared to in 
previous sections. 

In Figure 1181 we show S7(C iv) and S1(C) (or more ac- 
curately the carbon in the IGM) for the cw (top) and vzw 
(bottom) models. The cw model shows poor convergence, 
with more metals and C iv absorption being produced in 
higher resolution simulations. This is in stark contrast with 
the global star formation rate in these models, which shows 
excellent convergence at late times as demonstrated by SH03 
(higher-resolution simulations do resolve earlier star forma- 
tion). In the cw model, all galaxies roughly enrich the same 
volume independent of galaxy mass since they eject winds 
at the same speed. Hence higher resolution simulations that 
resolve star formation in smaller halos at earlier times dis- 
tribute metals much more widely. To obtain a convergent 
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Figure 18. The numerical resolution dependence of O(C) (as 
measured along our lines of sight) and f2(C iv) in the cw and vzw 
models in the 8, 16, and 32 /i -1 Mpc boxes with 256 3 gas and 
dark particles. The same data points described in Figure Il2l are 
displayed. f2(C) and H(Civ) converge quite well in momentum- 
driven wind models, unlike for the constant wind case. The 
C(Civ) measurement suggests slightly larger ionization correc- 
tions at lower redshift in larger boxes. 



IGM enrichment in a constant wind model requires resolv- 
ing every galaxy that has ever driven a wind, which is excep- 
tionally challenging computationally. Since the disagreement 
with observations is even greater in higher resolution runs, 
this model is unlikely to be successful in explaining IGM 
enrichment. 

The vzw model, on the other hand, shows much better 
resolution convergence over most of the redshift range in 
both n(Civ) and fi(C), with a minor exception at z > 5 
where the stochastic nature of biased early star formation 
may be producing differences between the realizations. 

The impact of a finite volume is also seen in the vzw 
model at z < 3. Here, the 32/i _1 Mpc volume produces a 
larger global Civ ionization correction, owing to larger and 
hotter structure being formed in the larger volume. The 
differences are small, but illustrate that large volumes are 
needed as well as high resolution in order to properly model 
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the IGM. Hence it is desirable to use the largest possible 
volume that is resolution-converged. 

Despite these minor issues, our overall conclusions 
should be robust as long as the simulations are capable of 
accurately resolving the properties of C IV absorbers. As Fig- 
ure ED shows, this is a good assumption for the CDD, with 
the 8 and 16/i _1 Mpc vzw volumes showing good convergence 
over most of the redshift range. The other momentum-driven 
winds follow similar trends, while cw is not converged (as 
with Q(C iv); not shown). The deviations in the vzw case are 
primarily associated with the 32/i -1 Mpc volume. It under- 
estimates weak absorbers at high-z because it cannot fully 
resolve all metal production, and it underestimates strong 
absorbers at low-z because strong C IV systems arise in small 
galactic halos that it cannot resolve. Hence our 16 /T^Mpc 
volume is optimal. Unfortunately, we cannot directly deter- 
mine if our 16 ft _1 Mpc volume is robust at z < 3, since we 
have only evolved our 8 h, _1 Mpc volume to z — 3 due to 
CPU time considerations. Some effects of multi-phase col- 
lapse in high density regions may therefore persist in this 
volume, which may explain the deficit present in all models 
compared to data at the largest column densities. 

The linewidths provide an even more stringent test of 
resolution, because they are keenly sensitive to small fluctu- 
ations in diffuse gas. Figure [50] shows the fe-parameter distri- 
bution for the vzw model in three redshift intervals. While 
the 8 and 16/i _1 Mpc agree very well, the w32n256vzw model 
clearly shows wider fc-parameters at all redshifts. If we split 
the width of the lines up into their thermal and spatial com- 
ponent as we did in H5.2.2I we find that the increase is mostly 
in b sp . At z > 4.5, we find w(Civ)= 12.7, 11.5, and 16.5 kpc 
for the 8, 16, and 32 h~ Mpc boxes respectively. Hence the 
line widths, as with the CDDs, are well converged up to 
16/i _1 Mpc but falter at lower resolution. 

Can we understand why the 16 h~ Mpc volumes are 
converged but 32 fa -1 Mpc ones are not? To do so, we es- 
timate the minimum spatial resolution needed to properly 
model IGM fluctuations at any given density. We begin 
with the minimum SPH resolution, which is approximately 
0.5 lbox/n P art x {p/p)~ 1/3 , or 0.5 d m is{p/ p)~ 1/3 where d mis 
is the mean interparticle spacing (i.e. a single resolution el- 
ement generally requires two SPH particles in any given di- 
rection). The physical minimum optically-thin absorber size 
depends on the Jeans length, the limit below which the ther- 
mal kinetic energy prevents the IGM from collapsing. The 
Jeans length is Tj = 975 yjTi/[{p/ p){\ + z) 3 ] kpc at a given 
density p, where T4 = T/10 4 . Equating the minimum SPH 
resolution to rj, we derive the following expression for the 
maximum interparticle spacing allowable in a simulation to 
resolve an absorber, 

d re8 = 170 (p/p)- 1/6 (4/(1 + z)y' i/2 VTl h^kpc. (7) 

Because of the inverse density dependence in d res , SPH 
simulations actually tend to resolve the low-density IGM 
fluctuations slightly better than the high-density ones, at a 
given temperature. However, particularly at high redshifts, 
low-density fluctuations tend to be colder (see Figure ITU . 
making them the dominant resolution constraint. The typi- 
cal value of 170 h~ 1 kpc above is uncomfortably close to our 
32 /i _1 Mpc volume's d m is = 125/i _1 kpc, hence this volume 
is starting to fail to resolve some absorbers, resulting in poor 
convergence. 



At low-z, ISimcoe et alJ i200rl) observed strong metal- 
line absorbers to have T ~ 10 4 K and p/p ~ 100 between 
z — 2.7 — > 1.8. These are similar to the physical parameters 
of our strongest C IV systems at these redshifts. At z — 2, we 
obtain d res — 51/i~ 1 kpc for these parameters, which would 
not quite be resolved even in our 16 ft -1 Mpc runs {d m %s = 
63/i _1 kpc). Hence the deficit in strong systems in our models 
may be due to lack of resolution. 

To summarize, the resolution convergence for the vzw 
model is quite good up to our 16h~ Mpc volume, which is 
precisely why we chose to quote predictions for this volume 
throughout the paper. Our 32/i _1 Mpc shows poorer conver- 
gence, but it also suggests that at lower redshifts (z < 3) 
the small 16/i _1 Mpc volume may slightly impact strong ab- 
sorber statistics. We have checked that our results for the 
vzw case are generally applicable to all the momentum- 
driven wind models. In contrast, the cw model shows poor 
resolution convergence, and the disagreements with obser- 
vations are significantly exacerbated as one moves to higher 
resolution. 



7 CONCLUSIONS 

We have explored various models of galactic outflows in- 
corporated into cosmological hydrodynamic simulations of 
structure formation, and compared them with observations 
of IGM enrichment at z > 1.5. With n o outflows, we show 
(in agreement with previous studies, e.g. lAguirre et all200ll) 
that a smattering of metals are distributed out to moder- 
ate overdensities through dynamical interactions, but can- 
not come close to enriching the diffuse IGM to observed 
values. Hence high-velocity outflows are required to enrich 
the IGM to observed values. In order to allow flexibility and 
perform controlled tests of various scenarios, we implement 
these outflows "by hand", tying the wind velocities (v w ) and 
mass loading factors (77) to galaxy properties through vari- 
ous parameterizations. 

In this paper we show tests of five outflow models. 
Most of these models, and all successful ones, are based on 
momentum-driven wind scenarios that are favored by obser- 
vations of local starburst galaxies. We consider a constant- 
wind (cw) model that produces strong, early outflows; a 
weak momentum-wind model (mw) that yields little IGM 
heating; an "observed" wind model (zw) obtained by taking 
the observed scalings of local outflows at face value (not 
a momentum-driven scenario); a strong momentum- wind 
model ( mzw) calibrated by local starburst observations of 
iMartinl d2005T) ; and a variable momentum wind model (vzw) , 
where we have incl uded a spread in wind speeds in a ccord 
with observations of lRupke. Veilleux fc Sanders! |2005). Our 
models are chosen to span an interesting range in outflow 
strengths, while broadly reproducing the global star forma- 
tion rate density evolution down to z — 1.5. 

Our first main result is that the outflow models favored 
by local starburst observations, when applied to galaxies 
forming hierarchically across cosmic time, are able to repro- 
duce a wide range of C IV absorption line observations dur- 
ing the epoch of peak cosmic star formation (z ~ 1.5 — 5.5). 
We further demonstrate that this agreement is non-trivial, 
requiring particular levels of metal and energy input into 
the IGM as a function of time. In particular, winds should 
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Figure 19. Column density distributions of C IV for the w8n256vzw, wl6n256vzw, and w32n256vzw simulations, compared to the BSR03 
points (black triangles). The 8 and 16 h~ 1 Mpc simulations appear resolution converged, while the 32 /i — 1 Mpc box underestimates weak 
absorbers at high-2 and possibly stronger absorbers at low-z. 
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Figure 20. The ft-parameters distributions in our three volumes for the vzw model. This demonstrates the w8n256vzw and wl6n256vzw 
absorbers are resolution-converged and fit the observations by BSR03, while the w32n256vzw absorbers are too wide. Using a Jeans 
length argument (see text), we show that the 32 h — x Mpc often has too poor resolution to properly resolve diffuse Civ absorbers. 



have high enough speeds to enrich the diffuse IGM at early 
times, but low enough speeds not to overheat the IGM. 
Winds should also have mass loading factors high enough 
to suppress early star formation, but low enough to produce 
enough early metals. Generically, momentum-driven wind 
models produce the desired scalings, because they heavily 
suppress early star formation but are able to distribute many 
of the metals formed into the IGM due to the high early 
mass loading factors in small galaxies. Yet, their low wind 
speeds from small galaxies do not overheat the IGM. It is 
therefore quite remarkable and non-trivial that our stronger 
momentum-driven wind scenarios (mzw and vzw), based on 



physical models calibrated by local observations of outflows, 
naturally reproduce a wide variety of observations. 

Our second key result is that the lack of evolution ob- 
served in the properties of C IV absorption such as the global 
mass density in Civ systems (f2(Clv) from z ~ 5.5 — * 1.5 
is not necessarily indicative of a lack of c ontinual enrich- 
ment, as has been argu ed previously (e.g. ISoneailal 1200 ll : 
IScannapieco et alj20 02). In all our models, most of the met- 
als are ejected into the IGM during the heydey of star forma- 
tion at z < 6, but our successful models also inject sufficient 
heat into the IGM to generate a drop in the C IV ionization 
fraction that roughly balances the increase in metallicity. 
The success of the mzw and vzw models in matching Q(C iv) 
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evolution renders unnecessary a putative population of early 
star formation that enriches a large fraction of the IGM at 

Comparing various models to f2(C iv) data showed that 
C IV is overproduced in the cw and zw models, primarily 
because early star formation is insufficiently suppressed and 
copious metals are produced. While the large heat input 
from the cw model partially counteracts metal overproduc- 
tion by lowering the C IV ionization fraction, it is insufficient 
to bring this model into agreement with data. Note that a 
co mparison to observ ed z ~ 6 galaxy luminosity functions 
bv lDave et alJ (120061) also indicated that star formation is 
not suppressed sufficiently early on in the cw model. One 
could envision lowering the wind speed in the cw model to 
reduce the heat input and lower fi(Civ), but this would 
exacerbate the overcooling problem at high redshifts. The 
highest- z observations of f2(C iv) tend to favor models that 
suppress early star formation, which occurs naturally in the 
momentum-wind driven scenarios (mzw, vzw) owing to high 
mass loading factors in small galaxies. 

The column density distributions of C iv absorbers 
present a test of the spatial distribution of metals. Mod- 
els with overly weak winds such as mw cannot pollute the 
low overdensity (p/p ~ few) IGM to reproduce the popula- 
tion of weak absorbers seen at z ~ 3 — 4. Once winds are 
strong enough (vzw, mzw, even cw and zw), the shape and 
(lack of) evolution of the CDD is roughly as observed. 

Observations of the C iv linewidth (6-parameter) distri- 
bution shows little evolution, which again is generally repro- 
duced in all our models. The linewidths provide constraints 
on the temperature of absorbing gas, which seem to accom- 
modate some but not too much feedback energy input, as 
provided by the vzw, mzw, and mw models. 

Pixel optical depth (POD) ratios provide a different way 
of characterizing Civ absorption, yet yield broadly similar 
conclusions. The Civ/Hi POD ratio tests the amount and 
spatial distribution of the metals, because H I absorption is 
correlated with underlying gas density. Comparisons with 
our models show that cw distributes metals too widely, mw 
not widely enough, and mzw and vzw lying in the interme- 
diate range that agrees with data. The C m/C iv POD ratio 
is a sensitive measure of temperature, and shows that the 
cw and possibly even mzw heat the IGM too much, while 
the vzw and mw models match well. 

In short, vzw and to a slightly lesser extent mzw are the 
only models that match all observations of C iv absorption 
we have considered, while reproducing a promising match to 
the global observed star formation history. These observa- 
tions provide sensitive tests of the metallicity and temper- 
ature structure of the IGM, making the agreement of these 
models an impressive achievement. The fact that these mod- 
els are consistent with local observations of outflows pro- 
vides an interesting connection between starbursts-driven 
outflows at all epochs. 

Our simulations also provide insights into the physical 
effects of outflows on the IGM. Both the vzw and mzw mod- 
els show relatively low volume filling factor of Z > 1O _3 Z0 
gas, increasing from ~ 2 — > 10% from z = 6 — > 2. Hence 
widespread enrichment of the IGM is not required to match 
observations, and in fact is disfavored as models that dis- 
tribute metals more widely tend to overproduce C iv absorp- 
tion. These results broadly agree with measurements of the 



quiescence of the IGM at high redshift that sug gest outflows 
cann ot have impacted a significant volume iRauch et alJ 

2005) . The metallicity-overdensity relation in the IGM is 
not a power-law as is sometimes assumed, and instead shows 
relatively fiat enrichment at high overdensities, dropping off 
steeply at lower overdensities. An interesting prediction of 
our models is that a significant fraction of C iv lines are colli- 
sionally ionized. However, this does not violate temperature 
constraints on the IGM as seen from &-parameter distribu- 
tions or Cm/Civ POD ratios. The median Civ absorber 
traces higher overdensity and hotter gas at lower redshifts, 
and the optical depth is only weakly correlated with under- 
lying gas density, u nlike the tight re lationship inferred for 
Lya absorbers (e.g. iDave et alJll999l) . This primarily owes 
to the large range of ionization fractions in C iv systems, 
from cooler photoionized absorbers to collisionally ionized 
systems reaching T ~ 3 x 10 K. The idea that C iv absorp- 
tion arises almost exclusively in diffuse, photoionized gas 
following a tight density-temperature relation similar to H I 
absorbers is not favored by our models. 

While we present five wind models in this paper, we 
have actually simulated over a dozen wind models, and de- 
cided to focus on these five because they are physically 
and/or observationally motivated, and best illustrate the re- 
lationship between outflow properties and IGM enrichment. 
Nevertheless, the range of outflow models investigated in 
this work is far from exhaustive, as one could envision more 
complex recipes that might still be consistent with obser- 
vations. Hence our favored models should be regarded as 
plausible, but not uniquely successful. 

Further uncertainties in the exact model that best de- 
scribes galactic outflows comes from a range of system- 
atic uncertainties we have not extensively considered in this 
work. One uncertainty is the impact of Helium reionization 
(investigated briefly in !j5. II . which can impact the com- 
parison of C iv absorption at early epochs. Unfortunately 
including this effect requires fairly precise knowledge of the 
topology and epoch of Helium reionization, which is poorly 
constrained at present. A second uncertainty is that our sim- 
ulations were run with a cosmology that is slightly different 
than that now favored by the 3rd-year WMAP release. In 
particular, the new values of ag (now 0.75 as opposed to 0.9 
as we used) and Q m (0.24 vs. 0.30) will signi ficantly lower 
the h alo collapse fraction at high redshift dSpergel et alJ 

2006) . This means that the mass loading factors will need 
to be reduced relative to equation |3] to match observations 
of high-redshift star formation rates. The net effect should 
be small: Given that observations dictate a certain global 
star formation rate, any wind model that matches it will 
produce a similar amount of metals. A third uncertainty in- 
volves the possible contribution of Population III stars to 
early IGM enrichment. Our models, while not requiring it, 
cannot rule out some contribution from early massive star 
formation that might have higher yields and larger super- 
nova energy input. More generally we have not considered 
the pos sibility of IMF vari ations, although recent work hints 
at this iFardal et al.l2006T) . A fourth uncertainty arises from 
possible contributions from AGN feedback in unbinding en- 
riched gas from galaxies, as bri ght quasars have been ob- 
served at z > 6 iFan et al.ll2006h . A fifth uncertainty is the 
effect of enriched stellar winds tecycling metals into the ISM, 
mostly from asymptotic giant branch stars, which will in- 
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crease the net carbon yield from stars. Despite these caveats, 
the fact that a plausible feedback model exists connecting 
local starbursts and high-z IGM enrichment represents sig- 
nificant progress towards quantifying the impact of outflows 
on galaxies and the IGM across cosmic time. 

Further constraints on outflows can be provided by ex- 
amining IGM metal evolution to z ~ 0, studying early IGM 
metal enrichment extending into the reionization epoch, 
matching the mass-metallicity evolution in galaxies, and re- 
producing the abundances patterns and gradients in galaxy 
clusters, among other things. We are currently investigat- 
ing all these aspects, especially the imp act of various out - 
flow models on galaxy populations (see iDave et alJ 12006^ . 
On the theoretical side, it is of course preferable to gener- 
ate and drive winds self-consistently out of galaxies rather 
than having to include them "by hand" , and progress is be - 
ing made towards that end (e.g. IScannapieco et all [2006). 
Observationally, we stress the need to obtain z > 4.5 ob- 
servations of C IV capable of resolving thermal linewidths in 
order to constrain the physical nature of early feedback and 
metallicity enrichment. In short, there remains much work 
to be done in order to obtain a fully self-consistent model of 
metal evolution and distribution in the universe. Our work 
represents a small but important step forward in this area. 



ACKNOWLEDGMENTS 

We thank A. Aguirre, K. Finlator, L. Hernquist, N. Katz, M. 
Pettini, E. Ryan- Weber, J. Schaye, and R. Simcoe for their 
helpful input. We thank V. Springel and L. Hernquist for 
the use of Gadget-2 before public release, and E. Scan- 
napieco for providing us with his database of C IV lines. 
The simulations were run on the Xeon Linux Supercluster 
at the National Center for Supercomputing Applications, 
and on our department's 100-processor Beowulf system at 
the University of Arizona. Support for this work, part of 
the Spitzer Space Telescope Theoretical Research Program, 
was provided by NASA through a contract issued by the 
Jet Propulsion Laboratory, California Institute of Technol- 
ogy under a contract with NASA. Support for this work 
was also provided by NASA through grant number HST- 
AR-10647.01 from the SPACE TELESCOPE SCIENCE IN- 
STITUTE, which is operated by AURA, Inc. under NASA 
contract NAS5-26555. 



REFERENCES 

Adelberger, K. L., Steidel, C. C, Shapley, A. E., & Pettini, 

M. 2003, ApJ, 584, 45 
Adelberger K. L., Shapley A. E., Steidel C. C, Pettini M., 

Erb D. K., & Reddy N. A. 2005, ApJ, 629, 636 
Aguirre, A., Hernquist, L., Schaye, J., Katz, N., Weinberg, 

D. H., & Gardner, J. 2001, ApJ, 561, 521 
Aguirre, A., Schaye, J., & Theuns, T. 2002, ApJ , 576, 1 
Aguirre, A., Schaye, J. Hernquist, L., Kay, S., Springel, V., 

& Theuns, T. 2005, ApJ, 620, L13 
Aracil, B., Petitjean, P., Picho, C, & Bergeron, J. 2004, 

A&A, 419, 811 
Balogh, M. L., Pearce, F. R., Bower, R. G., & Kay, S. T. 

2001, MNRAS, 326, 1228 



Becker, G. D., Sargent, W. L. W., Rauch, M., & Simcoe, 

R. A. 2006, ApJ, 640, 69 
Bertone, S. & White, S. D. M. 2006, MNRAS, 367, 247 
Boksenberg, A., Sargent, W. L. W., & Rauch, M. 2003, 

ASP Conference Proceedings, Vol. 297, 447, eds. Edited 

E. Perez, R.M.G. Delgado, & G. Tenorio-Tagle (BSR03) 
Bolton, J. S., Haehnelt, M. G., Viel, M., & Carswell, R. F. 

2006, MNRAS, 366, 1378 
Carswell, R. F., Webb, J. K., Baldwin, J. A., & Atwood, 

B. 1987, ApJ, 319, 709 

Cen, R., Nagamine, K., & Ostriker, J. P. 2005, ApJ, 635, 
86 

Cole, S., et al. 2001, MNRAS, 326, 255 

Croft, R. A. C, Weinberg, D. H., Katz, N., & Herquist, L. 

1998, ApJ, 495, 44 
Dave, R., Hernquist, L., Weinberg, D. H., & Katz, N. 1997, 

ApJ, 477, 21 

Dave, R., Hellsten, U., Hernquist, L., Katz, N., & Wein- 
berg, D. H. 1998, ApJ, 509, 661 
Dave, R., Hernquist, L., Katz, N., & Weinberg, D. H. 1999, 

ApJ, 511, 521 
Dave, R. et al. 2001, ApJ, 552, 473 
Dave, R. & Tripp, T. M. 2001, ApJ, 553, 528 
Dave, R., Finla tor, K., fc Oppenheimer, B. D. 2006, MN- 
RAS, in press, astro-ph/0511532 
Dekel, A. & Silk, J. 1986, ApJ, 303, 39 
Eisenstein, D. J. & Hu, W. 1999, ApJ, 511, 5 
Ellison, S. L., Pettini, M., Lewis, G. F., Songaila, A., Cowie, 

L. L. 1999, Ap&SS, 269, 201 
Ellison, S. L., Songaila, A., Schaye, J., & Pettini, M. 2000, 

AJ, 120, 1175 
Engelbracht, C. W. et al. 2006, ApJ, 642, L127 
Erb, D. K., Shapley, A. E., Pettini, M., Steidel, C. C, 
Red dy, N. A, k, A delberger, K. L. 2006, ApJ, accepted, 
|astro-ph/ 0602473 
Fan, X. et al. 2006, A.J. 132, 117 

Fardal, M. A., Katz, N., Weinberg, D. H., & Dav'e, R. 2006, 

MNRAS, submitted, astro-ph/0604534 
Ferrara, A., Pettini, M., & Shchekinov, Y. 2000, MNRAS, 

319, 539 

Fujita, A., Mac Low, M.-M., Ferrara, A., Meiksin, A. 2004, 
ApJ, 613, 159 

Frye, B. L., Tripp, T. M., Bowen, D. B., Jenkins, E. B., & 
Sembach, K. R. 2003, in "The IGM/Galaxy Connection: 
The Distribution of Baryons at z=0", ASSL Conference 
Proceedings Vol. 281, 231, eds. J.L. Rosenberg & M.E. 
Putman 

Gnedin, N. Y. & Ostriker, J. P. 1997, ApJ, 486, 581 

Haardt, F. & Madau, P. 1996, ApJ, 461, 20 

Haardt, F. & Madau, P. 2001, in "Clusters of galaxies and 
the high redshift universe observed in X-rays, Recent re- 
sults of XMM-Newton and Chandra" , XXXVIth Rencon- 
tres de Moriond , eds. D.M. Neumann & J.T.T. Van. 

Heckman, T. M., Lehnert, M. D., Strickland, D. K., & Ar- 
mus, L. 2000, ApJS, 129, 493 

Heckman, T. M. 2003, in Avila-Reese V., Firmani C, Frenk 

C. S., Allen C, eds, Rev. Mex. Astron. Assoc., Vol. 17, 
Galaxy Evolution: Theory & Observations, p. 47 

Hellsten, U., Hernquist, L., Katz, N., & Weinberg, D. H. 

1998, ApJ, 499, 172 
Hernquist, L., Katz, N., Weinberg, D. H., & Miralda- 

Escude, J. 1996, ApJ, 457, :51 



Simulations of IGM Enrichment 31 



Hui, L. & Gnedin, N. Y. 1997, MNRAS, 292, 27 
Hopkins, A. M. 2004, ApJ, 615, 209 

Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., 
Cole, S., Evrard, A. E., Couchman, H. M. P., & Yoshida, 
N. 2001, MNRAS, 321, 372 

Katz, N., Weinberg, D. H., & Hernquist, L. 1996, ApJS, 
105, 19 

Kennicutt, R. C. 1998, ApJ, 498, 541 

Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, 

MNRAS, 363, 2 
Kirkman, D., Tytler, D., Suzuki, N., Melis, C, Hollywood, 

S., James, K., So, G., Lubin, D., Jena, T.,; Norman, M. 

L., & Paschos, P. 2005, MNRAS, 360, 1373 
Madau, P., Ferguson, H. C, Dickinson, M. E., Giavalisco, 

M., Steidel, C. C, & Fruchter, A. 1996, MNRAS, 283, 

1388 

Mailer, A. H., & Bullock, J. S. 2004, MNRAS, 355, 694 
Martin, C. L. 1999, ApJ, 513, 156 
Martin, C. L. 2005, ApJ, 621, 227 
McKee, C. F. & Ostriker, J. P. 1977, ApJ, 218, 148 
Mo, H. J. & Miralda-Escude, J. 1996, ApJ, 469, 589 
Murray, N., Quatert, E., & Thompson, T. A. 2005, ApJ, 
618, 569 

Olbers, H. W. M. 1826 

Page et al. 2006, ApJ, submitted, | |astro-ph/ 0603450 
Pettini M., Shapley A. E., Steidel C. C, Cuby J.-G., Dick- 
inson M., Moorwood A. F. M., Adelberger K. L., & Gi- 
avalisco M. 2001, ApJ, 554, 981 
Pettini, M., Madau, P., Bolte, M., Prochaska, J.X., Ellison, 

S.L., & Fan, X. 2003, ApJ, 594, 695 
Porciani, C. & Madau, P. 2005, ApJ, 625, L43 
Prochaska, J. X., Gawiser, E., Wolfe, A. M., Castro, S., & 

Djorgovski, S. G. 2003, ApJ, 595, L9 
Press, W. H., Rybicki, G. B., & Schneider, D. P. 1993, ApJ, 
414, 64 

Rauch, M., Haehnelt, M. G., & Steinmetz, M. 1997, ApJ, 
481, 601 

Rauch, M., Sargent, W. L. W., & Barlow, T. A. 2001, ApJ, 
554, 823 

Rauch, M., Becker, G. D., Viel, M., Sargent, W. L. W., 
Smette, A., Simcoe, R. A., Barlow, T. A., Haehnelt, M. 

G. 2005, ApJ, 632, 58 

Rupke, D. S., Veilleux, S., & Sanders, D. B. 2005, ApJS, 
160, 115 

Ryan- Weber, E. V., Pettini, M., fc Madau, P. 2006, MN- 
RAS, accepted, astro-ph/0607029 

Scannapieco, E., Ferrara, A., & Madau, P. 2002, ApJ, 574, 
590 

Scannapieco, E., Pichon, C, Aracil, B., Petitjean, P., 
Thacker, R.J., Pogosyan, D., Bergeron, J., & Couchman, 

H. M.P. 2006, MNRAS, 365, 615 
Schaerer, D. 2003, A&A, 397, 527 

Schaye, J., Theuns, T., Rauch, M., Efstathiou, G., & Sar- 
gent, W. L. W. 2000, MNRAS, 318, 817 

Schaye, J., Aguirre, A., Kim, T.-S., Theuns, T., Rauch, M., 
& Sargent, W.L.W. 2003, ApJ, 596, 768 (S03) 

Shapley, A. E., Steidel, C. C, Pettini, M., & Adelberger, 
K. L. 2003, ApJ, 588, 65 

Simcoe, R.A., Sargent, W.L.W., Rauch, M., & Becker, G. 
2006, ApJ, 637, 648 

Simcoe, R. A. 2006, ApJ, submitted, astro-ph/0605710 

Sokasian, A., Abel, T., & Hernquist, L. 2003, MNRAS, 340, 



473 

Songaila, A., Cowie, L.L. 1996, AJ, 112, 335 
Songaila, A. 1998, AJ, 115, 2184 
Songaila, A. 2001, ApJ, 561, L153 

Songaila, A. 2005, AJ, 130, 1996 

Spergel et al. 2006, ApJ, submitted, astro-ph/0603449 
Springel, V. & Hernquist, L. 2002, MNRAS, 333, 649 
Springel, V. & Hernquist, L. 2003, MNRAS, 339, 289 
Springel, V. & Hernquist, L. 2003, MNRAS, 339, 312 
(SH03) 

Springel, V. 2005, MNRAS, 364, 1105 
Theuns, T., Schaye, J., & Haehnelt, M. G. 2000, MNRAS, 
315, 600 

Theuns, T., Bernardi, M., Frieman, J., Hewett, P., Schaye, 

J., Sheth, R. K., & Subbarao, M. 2002, ApJ, 574, Llll 
Theuns, T., Viel, M., Kay, S., Schaye, J., Carswell, R.F., 

& Tzanavaris, P. 2002, ApJ, 578, L5 
Tremonti, C. A., et al. 2004, ApJ, 613, 898 
Sutherland, R. S. & Dopita, M. A. 1993, ApJS, 88, 253 
Woosley, S. E. & Weaver, T. A. 1995, ApJS, 101, 181 
Zheng, W., et al. 2004, ApJ, 605, 631 



